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Abstract 


In this study I use partial frequency redistribution (PRD) functions to examine their effects on 
spectral line formation in spherically symmetric and expanding atmospheres of cool giant and 
supergiant stars. Primary aim of this investigation is to bring out the differences between the 
emergent spectral line profiles resulting from PRD and complete redistribution (CRD) under 
the influence of various physical parameters characterizing the atmospheres of cool giants and 
supergiants. The appreciation of this aspect will be important for quantitative analysis of 
stellar spectra and for computing model atmospheres of such stars. 

In a scattering process, both the direction and the frequency of a photon may change. These 
changes are described by partial frequency redistribution functions. There are following four 
categories of redistribution functions: Case I, zero line width, denoted by Rj. It does not 
apply to any real line; demonstrates the effects of Doppler redistribution alone to an observer 
in the laboratory frame. Case II, radiation .damping in the upper state and coherence in 
the atom*s rest frame, denoted by Rij. It applies to resonance lines in low density media. 
Case in, complete redistribution in the atom*s frame, denoted by J?///. Case IV, subordinate 
line redistribution between two broadened'states, denoted by Jf2v> It applies to non-coherent 
subordinate line scattering. Complete redistribution is a limiting approximation which implies 
that there is a complete reshufiUng of atoms in their excited state in such a way that there is 
no correlation between the frequencies of the incoming and the scattered photons. 

The assumption of CRD has been widely used in earlier works on line transfer because it 
not only simplifies the numerical solution of the transfer equation but also provides a good 
approximation to reality in those media which are dense enough to support high collision rates. 
However, in the extended and tenuous atmospheres of cool giants and supergiants where low 
densities (and hence low collision rates) prevail, CRD is not expected to yield accurate emergent 
spectral line profiles. Moreover, in the wings of strong resonance lines, the scattering is nearly 
coherent. This leads to a deviation of the line profiles from those calculated in accordance with 
the assumption of CRD. Therefore, a natural recourse to the application of partial frequency 
redistribution can lead to more accurate line profiles. 

Over the last two decades, considerable progress has been made in using partial frequency 
redistribution functions to study the spectral line formation in idealized stellar atmospheres. 
Most of these studies are limited to the assumption of plane-parallel geometry and/or the 
absence of velocity fields. These assumptions are unrealistic because the real atmospheres of 
cool giants and supergiants expand and are geometrically extended. There are only few papers 
in literature which have taken into account both the spherical geometry and the expansion 
effects to study PRD effects on emergent line profiles in idealized atmospheres. These studies 
employ either Rj or J?//. Rm has usually been represented by CRD. Some studies have been 
done to model the spectral lines formed in the expanding chromospheres of red giants using 
jR/j. So far, the combined effects of sphericity and velocity fields on the differences between the 
solutions resulting from Ry and the CRD have not been explored. As an important feature 
of the present work I make a detailed comparative study of the PRD effects of Rn, 



Riii and Rv on the spectral line formation in the atmospheres of cool giants and 
supergiants taking into account both the sphericity and the expansion effects. 

hi order to achieve the above mentioned goals, I adopt a two-pronged approach. In the first 
approach, I consider parameterized stellar atmospheres to explore the influence of different 
physical parameters on the differences between the solutions due to PRD and CRD. These 
include the atmospheric extent, expansion velocity and the thermalization parameter (which 
gives the probability of coUisional destruction of photons). The solutions considering idealized 
models serve as a preliminary step to solutions of actual stellar atmospheric problems by 
providing an insight into the basic physical phenomena. 

hi the second approach, I adopt computed model atmospheres (temperature and electron 
number density distributions) of red giants and supergiants. Prom these models, the line optical 
depth and the thermalization parameter is computed at each depth point in the atmosphere 
to study the frequency redistribution effects on some actually observed spectral lines from the 
stars of the same type. This is an important study as it brings out the model dependence of the 
differences between the partial frequency redistribution and complete redistribution solutions. 

For all cases, I solve the equation of transfer in the frame work of discrete space theory of 
radiative transfer using Peraiah-Grant method. 

In the first parametric study I examine the effects of angle-averaged Rv and Rn on emergent 
flux profiles in static spherically symmetric stellar envelopes. I write the equation of transfer 
for spherical symmetry in the rest frame for a two-level atom. Initially, to bring out the 
effects of geometrical extension alone, I neglect velocity fields. I assume a constant value of the 
thermalization parameter and set Planck’s function equal to unity throughout the atmosphere. 

I neglect the background continuum. In this study I show that the effect of spherical extension 
is much more on the emergent flux profiles resulting from Rn than on those resulting from Ry 
when the damping parameters of the lower and upper levels of an atom are comparable. The 
effect of Rv (as compared to CRD) on emergent flux profiles increases by making the lower 
level sharper than the upper level. This effect is further enhanced by increasing the spherical 
extension. 

In the second parametric study I go a step further to include the effects of expansion 
velocity in order to explore the combined effects of sphericity and expansion on the emergent 
flux profiles resulting from Rn and Ry, as compared to CRD. I consider an atmosphere having 
geometrical extension and small velocity of expansion which is characteristic of the atmospheres 
of cool giants and supergiants. I assume that there is no background continuum; set Planck’s 
function equal to unity and choose a constant value of the thermalization parameter at all 
depths in the atmosphere. In this study I show that there are substantial differences between 
the emergent flux profiles resulting from PRD {Rn and Ry) and CRD. 

In the third parametric study I investigate the PRD effects of JRj/j alone on the emergent 
line profiles in an expanding spherically symmetric stellar atmosphere. This is because, the 
limiting case of complete redistribution is usually taken to be an adequate representation of 
Riii in radiative transfer problems to study the spectral line formation in stellar atmospheres. 



However, the combined effect of sphericity and macroscopic velocity cannot be estimated a 
priori. I consider only low velocity regime which is appropriate for the atmospheres of cool 
giants and superpants. I treat this problem in the rest frame formalism. In this work I 
established that the differences between the emergent mean intensity profiles resulting from 
CRD and Rm are much smaller in expanding spherically symmetric stellar atmospheres than 
reported previously for plane-parallel atmospheres. 

As a next important step in this investigation, I adopt computed model atmospheres of red 
giants and supergiants to explore the PRD effects (as compared to CRD) on some observed lines 
(e.g. Si II). I use the temperature structure and the electron k hydrogen number densities from 
the given models to compute the line optical depth and the thermalization parameter at each 
depth point in the atmosphere. Model values of the radii and the atmospheric extensions are 
used. I solve the line transfer equation using Ru and CRD to compute the emergent spectral 
line profiles. This study shows that the differences between the results due to PRD and CRD 
are highly model dependent. The results of non-chromospheric models show a departure of 
«10%; where as the results of chromospheric,model of a-Orionis show enormous departures 
of PRD results from those of CRD. 

Finally, I summarize the results of this study along with its achievements and limitations. 
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Chapter 1 
Introduction 


The title Atmospheres of Cool Stars is very broad. No single study can cover all its aspects. 
This investigation focuses on one of the most important aspects of the subject : effects of 
partial frequency redistribution functions on the spectral line formation. Spectral 
lines carry a wealth of information concerning, the medium in which they form. The line 
spectrum can be used to obtain accurate empirical models of stellar atmospheres, complete with 
such quantities as temperature, gas pressure, chemical composition, gravitational acceleration, 
velocity fields. With the advent of observations from balloons, rockets and satellites, great 
improvements have been made in the quantity and quality of stellar spectroscopic data. The 
availability of highly accurate line profiles has enhanced the importance of partial frequency 
redistribution in the process of spectral line formation. 


1.1 Motivation 


In the process of scattering in spectral line formation, an atom is excited from one bound level 
to another by absorbing a photon, and then decays back radiatively to the original state, with 
the emission of a photon. Very often, in the works on line transfer, it is assumed that the 
photons are completely redistributed over the line profile. This assumption has been widely 
used because, it not only simplifies the numerical solution of the transfer equation but also 
provides a good approximation to reality in those media which are dense enough to support high 
collision rates. However, in the extended and tenuous atmospheres of cool giant and 
supergiant stars where low densities (and hence low collision rates) prevail, CRD is 
not expected to yield accurate emergent spectral line profiles. Thus, it is necessary to 
consider the redistribution of photons by employing partial frequency redistribution functions, 
which describe the scattering process precisely. Moreover, in the wings of strong resonance 
lines, the scattering is nearly coherent. This leads to a deviation of the line profiles from those 
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ealculated in accordaace with the assumption of CRD. Therefore, a natur 2 tl recourse to the 
application of partial frequency redistribution (PRD) can lead to more accurate line profiles. 


1.2 Brief definitions of partial frequency redistribution 
functions 


In a scattering process, both the direction and the frequency of a photon may change. These 
changes are described by partial frequency redistribution functions. There are following four 
categories of redistribution functions: Case I, zero line widths denoted by Rj. It does not 
apply to any real line; demonstrates the effects of Doppler redistribution alone to an observer 
in the laboratory frame. Case II, radiation damping in the upper state and coherence in 
the atom*s rest frame^ denoted by Rjj. It applies to resonance lines in low density media. 
Case III, complete redistribution in the atom*S:frame, denoted by Case IV, subordinate 
line redistribution between two broadened states^ denoted by Rv- It applies to non-coherent 
subordinate line scattering. Complete redistribution is a limiting approximation which implies 
that there is a complete reshuffling of atoms in their excited state in such a way that there is 
no correlation between the frequencies of the incoming and the scattered photons. 


1.3 Objective 


This study investigates the effects of partial frequency redistribution (PRD) functions Rj/, 
Riu and Rv on spectral line formation in spherically symmetric and expanding atmospheres 
of cool giant and supergiant stars. Primary objective of this investigation is to bring out 
the differences between the emergent spectral line profiles resulting from PRD and complete 
redistribution (CRD) under the influence of various physical parameters characterizing the 
atmospheres of cool giants and supergiants. As an important feature of this work a 
detailed comparative study is made of the PRD eflfects of Rj/, Rjjj and Ry on the 
spectral line formation in the atmospheres of cool giants and supergiants taking 
into account both the sphericity and the expansion effects. The appreciation of this 
aspect will be important for quantitative analysis of stellar spectra and for computing model 
atmospheres of such stars. 


1.4 Review of the related work 


Five decades have passed since Heney [27] discussed about the frequency redistribution of 
radiation including both the natural broadening and the Doppler motions of the radiating 
atom. Other early discussions on this subject are found in the works of Holstein [28], Unno 
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[69], Sobolev [67] and Warwick [76]. It was Hummer [31] who derived the redistribution 
functions for two levd atomic model. He defined the term redistribution function n*] u, n) 
as the probability that a photon having initial frequency i/ and direction n' will be absorbed, 
leading to the re-emission in the line of a (i/, n) photon. He then defined the following four 
redistribution functions : 

1. For the cases of zero line width in the atom’s frame; denoted by Rj. 

2. For radiative damping of the upper level only; denoted by Rn. 

3. For radiative and collisional damping of the upper level only; denoted by Rni- 

4. For radiative damping damping of both levels; denoted by Rjy. 

Hummer derived angle-independent redistribution functions for the cases of isotropic and dipole 
phase functions. He zdso discussed the symmetry properties of redistribution functions. In 
subsequent radiative transfer calculations, isotropic angle-independent redistribution functions 
have foimd mavimiiTn use. At this point it ig worth mentioning that, Heinzel [23] pointed 
out that Hummer’s [31] Rjv function is erroneous because it is derived firom an incorrect 
expression for scattering in the atomic frame. In the same paper Heinzel [23] derived a new 
redistribution function Ry for resonance scattering of subordinate lines assuming that both 
levels are radiatively broadened. Later, Heinzel and Hubeny [24] generalized this redistribution 
function to include collisional broadening of both the upper and lower levels of the atom. 
Jefferies and White [33] examined the validity of the CRD approximation and found that 
in the line core the photons are completely redistributed whereas in the wings scattering is 
nearly coherent. On this basis they proposed an approximation known as partially coherent 
scattering. Kneer [35] pointed out that this approximation is neither normalized nor symmetric, 
and he proposed a modified form of the partially coherent scattering redistribution function. 
Subsequently, Basri [9] found that partially coherent scattering cannot account for Doppler 
diffusion in the line wings, which can be important for low density atmospheres. Therefore, 
he argued that it is required to use exact forms of redistribution functions in order 
to achieve accurate radiative transfer solutions. Ayres [5] proposed a revised form of 
partially coherent scattering which includes the role of Doppler diffusion. Jefferies and White 
also called attention to the lowering of the line profile in the inner wing which occurs as the 
scattering in the wing is taken to be nearly coherent. They noted that this lowering can be 
caused in isothermal atmospheres by taking scattering to be coherent in the line wing, or in 
a chromosphere if the whole line fine is formed in CRD. Therefore, a proper understanding 
and treatment of spectral line formation in the line wings is essential for determining the 
temperature structure in a stellar atmosphere. 

Hummer [32] presented the first detailed solutions of the radiative transfer equation including 
PRD to go beyond the schematic partially coherent scattering approximation of Jefferies and 
White [33]. He used planar geometry, constant properties of the medium, two-level atomic 
model and angle-averaged Rj and R/j for isotropic phase function. He showed that the emer¬ 
gent intensity due to R/j in the line wings can be substzmtially lower than the CRD solutions. 
Vardavas [71, 72, 73] studied the effects of angle-averaged and angle-dependent R/, R/j and 
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Ru on spectral line formation in planar geometry for static and moving media. Peraiah [51] 
described the numerical solution of the transfer equation in spherical geometry bzused on the 
discrete-space theory of radiative transfer. He presented a computer code and the solutions of 
transfer equation in the rest frame using PRD functions Ri and Rn. Peraiah [52] used Rn 
to investigate its effects on spectral line formation in spherically S3nmmetric and expanding 
stellar atmospheres. Rangaxajan et al. [60] studied the effects of stimulated emission on rar 
diative transfer with PRD functions Ru and Rm in plane-parallel geometry. The combined 
effect of Rji and non-coherent electron scattering on polarized resonance line transfer in plane- 
parallel geometry was examined by Nagendra et al. [46]. Mohan Rao et al. [44] employed 
Domhe-Hubeny redistribution matrix and used a combination of Rji and Rni to study the 
polarized line transfer with coUisional redistribution in a plane-parallel atmosphere. Nagendra 
[47] also used Domke-Hubeny redistribution matrix to study PRD effects on resonance line 
polarization in spherical atmospheres. CoUisional redistribution effects on line polarization in 
spherical atmospheres were investigated by Nagendra [48]. Mohan Rao et al. [45] explored the 
effects of Ru, Rui and Ry on source functions {Sl) for plane-parallel and isothermal steUar 
atmospheres and showed, that, for a purely scattering medium with frequency independent 
incident radiation, the frequency dependent source function Si(Rv) lies below Sl{Ruj) but 
above Sl{Ru) in the line wings. Singh [65] examined the combined effect of sphericity and 
expansion velocity on the differences between the emergent mean intensity profiles resulting 
from Rui and CRD. 

Vernazza [74] showed that the flux and shape of the solar La line could be matched by a 
parameterized PRD calculation in which the core is computed with CRD assumption and the 
scattering in the wings was taken to be 93% coherent and 7% CRD. Omont, Smith and Cooper 
[49] made quantum mechanical calculations to derive the redistribution function for resonance 
lines. These results led Milkey and Mihalas [40] to propose a PRD computational scheme in 
which they solved the coupled equations of radiative transfer and statistical equilibrium using 
the Vernazza, Avrett, and Loeser (VAL) [75] solar model. They arrived at three important 
results: 1) A fuUy self-consistent PRD calculation using depth-dependent coherence fractions 
computed from the broadening rates is feasible; 2) The computed line wings lie between the 
CRD and pure coherent scattering limits and are similar to the 93% coherency calculations 
of Vernazza [74]; and 3) A proper inclusion of PRD effects in the wings of La is needed to 
compute electron density properly. In a subsequent paper Milkey and Mihalas [41] were able 
to match the shape of the observed La profile by including a depth-dependent line profile 
function. Roussel-Dupre [62] pointed out a major problems in the method of Milkey and 
Mihalas [40, 41]: the Omont et al. [49] redistribution function used by them is not valid for 
the La wings. Next, Milkey and Mihalas [42] applied their PRD method to the resonance 
line of Mgll treating each transition as a two-level atom. They showed that the line source 
function behaves like pure coherent scattering case, while in the line core it is independent of 
frequency and slightly above the CRD case. Milkey, Ayres, and Shine [39] examined the effect 
of gravity on the shapes of Mgll lines. They found that r^guiTsoco) increases systematically 
with decreasing gravity. Ayres, Linsky and Shine [6] and Ayres [4] used this result in trying 
to explain the Wilson-Bappu effect. 
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Shine, Milkey and Mihalaa [64] applied cross-redistribution formulation to a five-level, five- 
transition Call ion. They assumed HSRA model and several microturbulent velocity distri¬ 
butions. Their results provided strong evidence that the PRD approach includes much of the 
essential physics of the scattering processes, and in principle, PRD effects can account for many 
discrepancies between the line profiles calculated assuming CRD and the observations. Shine, 
Milkey and Mihalas [64] also considered the gravity dependence of the Call lines with results 
in qualitative agreement with the results of Milkey, Ayres, and Shine [39]. Milkey, Shine and 
Mihalas [43] then included the full angle and frequency dependence of redistribution during 
the scattering processes in an observer’s frame formulation. They foimd that angle-dependent 
effects are negligible for the Call line in a homogeneous solar chromospheric model, but these 
effects coidd be important for an inhomogeneous medium like the real solar chromosphere. 

Parallel with solar modeling, the development of chromospheric models for late-type stzurs based 
on the matching of observed and computed PRD profiles has taken place. Ayres and Linsky 
[7] constructed first such model to match Call K line of a Boo (Arcturus, K2 III). They could 
obtain a good agreement to the shape of the iqher wing and flux of the profile. In subsequent 
papers various authors who extended this work to derive atmospheric models for giant and 
supergiant cool stars are: Kelch et al [34] fitted K line of Call and h & k lines of Mgll using 
plane-parallel geometry, static medium for a Aur (G5 III), (3 Gem (KO III) and a Tau (K5 III); 
Baliunas et al [8] using plane-parallel geometry, static medium fitted the same line as Kelcl 
et al. [34] for A And (G8 III-IV) and a Aur (G5 III); Eriksson et al. [15] used plane-parallel 
geometry, static medium and fitted Mgll K line for (3 Cet (G9.5 III) star; Basri et al. [10] 
fitted Call K, Mgll K lines using plane-parallel static medium for /? Dra (G2 Il-Ib), c Gem 
(G8 Ib) and a Ori (M2 lab) stars; Basri et al. [10] fitted the same lines using plane-parallel 
geometry and circulation pattern for ^ Dra (G2 Il-Ib); Avrett and Johnson [3] used plane- 
parallel geometry, static medium to fit Mgll K line for N-type Carbon stars; Drake [14] used 
spherical geometry with winds to fit Mgll K line for a Boo (K2 III) star. In this list, except 
Basri [9], all other papers have neglected velocity fields and atmospheric extension. 

As appears evident from the review presented above, considerable progress has been made in 
using partial frequency redistribution functions to study the spectral line formation in idealized 
stellar atmospheres. But, most of these studies are limited to the assumption of plane-parallel 
geometry and/or the absence of velocity fields. These assumptions are unrealistic because 
the real atmospheres of cool giants and supergiants expand and are geometrically 
extended. Moreover, these studies employ either Rj or J7jj. Rm has usually been represented 
by CRD. So far, the combined effects of sphericity and velocity fields on the differences between 
the solutions resulting from Rji,Rv and the CRD have not been explored in a comparative 
way. It is this gap which this study tries to fill up. 


5 



1.5 Layout of the presentation 


In order to achieve the above mentioned goals, we adopt a two-pronged approach. In the 
first approach, we consider parameterized stellar atmospheres to explore the influence of 
different physical parameters on the differences between the solutions due to PRD and CRD. 
These include the atmospheric extent, expansion velocity and the thermalization parameter 
(which gives the probability of coUisional destruction of photons) . The solutions considering 
idealized models serve as a preliminary step to solutions of actual stellar atmospheric problems 
by providing an insight into the basic physical phenomena. 

In the second approach, we adopt computed model atmospheres (temperature and electron 
number density distributions) of red gieints and supergiants. From these models, the line optical 
depth and the thermalization parameter is computed at each depth point in the atmosphere 
to study the frequency redistribution effects on some actually observed spectral lines from the 
stars of the same type. This is an important study as it brings out the model dependence of the 
differences between the partial frequency redistribution and complete redistribution solutions. 

The work is presented in seven chapters. . 

• In chapter one, the need to study the partial frequency redistribution effects on the 
spectral line formation in the atmospheres of cool stars is discussed. Against this back¬ 
ground, relevant literature on the use of various partial frequency redistribution functions 
to study the spectral line formation in stellar atmospheres is reviewed. The chapter is 
concluded with the layout of the presentation. 

• In chapter two, the solution of equation of radiative transfer in the frame-work of the 
discrete space theory of raxiiative transfer using Peraiah-Grant formalism is described. 

• In chapter three, the effects of angle-averaged Rv and Rn on emergent flux profiles in 
static spherically symmetric stellar envelopes is examined. The equation of transfer for 
spherical symmetry in the rest frame is solved for a two-level atom. Initially, to bring 
out the effects of geometrical extension alone, velocity fields are neglected. A constant 
value of the thermalization parameter is assumed and the Planck’s function is set equal 
to unity throughout the atmosphere. The background continuum is neglected. This 
study establishes that sphericity enhances PRD effects. It is also shown that the effect 
of spherical extension is much more on the emergent flux profiles resulting from Rjj than 
on those resulting from Ry when the damping parameters of the lower and upper levels 
of an atom are comparable. The effect of Ry (as compared to CRD) on emergent flux 
profiles increases by making the lower level sharper than the upper level. This effect is 
further enhanced by increasing the spherical extension. 

• In chapter four, a step further is taken to include the effects of expansion velocity in 
order to explore the combined effects of sphericity and expansion on the emergent flux 
profiles resulting from Ru and Rv, as compared to CRD. An atmosphere having geomet¬ 
rical extension and small velocity of expansion is considered which is characteristic of the 


6 



atmospheres of cool giants and supergiants. It is assumed that there is no backgrotmd 
continuiim; Planck’s function is set equal to unity and a constant value of the thermal- 
ization parameter is chosen at all depths in the atmosphere. This study shows that there 
are substantial differences between the emergent flux profiles resulting from PRD {Rn 
and Rv) and CRD. 

• In chapter five, the PRD effects of i?/jj alone on the emergent line profiles in an expand¬ 
ing spherically symmetric stellar atmosphere are explored. This is because, the limiting 
case of complete redistribution is usually taken to be an adequate representation of R[jj 
itt radiative transfer problems to study the spectral line formation in stellar atmospheres. 
However, the combined effect of sphericity and macroscopic velocity cannot be estimated 
a priori. Only low velocity regime is considered which is appropriate for the atmospheres 
of cool giants and supergiants. This problem is treated in the rest frame formalism. In this 
work it is established that the differences between the emergent mean intensity profiles 
resulting from CRD and Rm are much smaller in expanding spherically s)mimetric stellar 
atmospheres than reported previously for plane-parallel atmospheres. 

• In chapter six, ais a next important step in this investigation, computed model atmo¬ 
spheres of red giants and supergiants are adopted to explore the PRD effects (as compared 
to CRD) on some observed lines (e.g. Si II). The temperature structure and the electron 
number density from the given models is used to compute the line optical depth and the 
thermalization parameter at each depth point in the atmosphere. Realistic values of the 
radii and the atmospheric extensions of cool giants are used. The line transfer equation 
is solved using Rn and CRD to compute the emergent spectral line profiles. This study 
aims at examining the model dependence of the differences between the results of PRD 
and CRD. 

• Finally, in chapter seven the highlights of the present work are summarized and possible 
directions for further study using partial frequency redistribution functions are indicated. 
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Chapter 2 

Solution of the line transfer equation 
using Discrete Space Theory 


2.1 Introduction 


The equation of radiative transfer, which describes the transport of radiation through a 
medium, is one of the essential equations in the stellar atmospheric research. This equa¬ 
tion incorporates the physical properties of the medium and its solution provides insight into 
the behavior of the internal and emergent radiation field. It is an. integro-differential equation 
and this feature increases the complexity of obtaining the solution depending upon the physics 
of the problem. It is not possible to solve this equation analytically for numerous astrophysical 
problems and so there exist several numerical methods to obtain the solution of the transfer 
equation. Amongst the notable methods, and the one used in this investigation is based on the 
Discrete Space Theory of Radiative Transfer. This method has two-fold advantages: 

1. The solution obtained can be subjected to rigorous numerical analysis to enable one 
to predict in advance the most economic finite difference mesh for a problem. 

2. The abstract mathematical structure is capable of generalization to other geometries. 

This chapter describes the basics of the Discrete Space Theory and the solution 
of line transfer equation based on this theory. 


2.2 Discrete Space Theory of Radiative Transfer 


The earliest reference to the Principle of Invariance, which is a special case of the Interaction 
Principle, is found in a paper by Stokes [68] in which he obtained difference equations for 
reflection by a pile of identical glass plates. However, it was Ambartzumian [2] who introduced 
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the Principle of Invaxiance into the radiative transfer theory. Chandrasekhar [13] extensively 
used these principles to calculate the reflection and transmission functions. The classical 
statement of the Principle of Invariance can be stated as: The law oj diffuse reflection of 
an infinitely homogeneous plane-parallel scattering medium is invariant with respect to the 
addition or subtraction of layers of arbitrary finite optical thickness to or from the medium. 
The Principle of Invariance is a consequence of the Interaction Principle and is essentially 
a statement of energy conservation in any medium. Redhefer [61] and Preisendorfer [59] wrote 
the Interaction Principle without the internal source terms. Grant and Hunt [19, 20] included 
the internal source terms and developed a numerical method to solve the transfer equation 
in inhomogeneous plane-parallel media. Grant and Peraiah [21] developed the solution of 
radiative transfer in spherical systems in the frame work of the Discrete Space Theory. This 
method has been applied to several case of extended and expanding stellar envelopes. 

Peraiah and Grant [57] presented a numerical method to obtain diffuse radiation in the ex¬ 
tended atmospheres of supergiant stars. This method is an extension of the discrete space 
theory of Grant and Hunt [19, 20] and has b^n widely used to study the spectral line for¬ 
mation in moving and extended atmospheres of supergiant stars. It has the following distinct 
merits: 


• It is a highly accurate method for computing the internal as well as emergent radiation 
in the presence or absence of the internal radiation sources. 

• It offers great flexibility to include a large number of physical processes that occur in 
stellar atmospheres. 

• It is numerically stable and gives physically meaningful solutions. 


In general the following steps are taken to obtain the solution 

1. The medium is divided into a number of cells whose thickness is less than or equal to 
the critical thickness (rcrit). The critical thickness is determined on the basis of the 
physical characteristics of the medium. Vcrit ensures stability and uniqueness of the 
solution, 

2. The integration of the transfer equation is performed on the cell which is a two- 
dimensional grid bounded by [r„,rn 4 i] x [/^j-i/aiA^i+i/a]) where 

J 

= S Cfc, j = 

fc=i 

where c* are the weights of Gauss Legendre quadrature formula. 

3. These discrete equations are compared with the canonical equations of the interaction 
principle smd to obtain the transmission and reflection operators of the cell. 

4. Lastly, all the cells are combined by star algorithm and obtain the radiation field. 
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The medium can be divided into shells whose thicknesses are larger thsm Tcrit but integration is 
done only on cell and star algorithm is used to obtain the transmission and reflection operators 
of the shells by adding the cells (star product). 

Further, a closer look at the Interaction Principle^ the Star Product and the Calculation of 
radiation field at interned points is presented, following Peraiah [56]. 


2.2.1 Interaction Principle 

The interaction principle relates the incident and emergent radiation field from a medium of 
given optical thickness. It permits us to set up exact difference equations for both reflection 
and transmission coefficients for radiation field. These coefficients are exact and can be used 
in the calculation of radiation field in the finite layers. 

Figure 2.1 shows a shell of optical thickness with incident and emergent intensities. It is 
assumed that specific intensities and U~^x incident at the boundaries n and n-f-l 
respectively of the shell with optical thickness r. The symbols with signs + and - represent 
specific intensities of the rays travelling in opposite directions. If /x represents the cosine of 
the angle made by a ray relative to the common normal (in the spherical case, /i is the angle 
made by the ray direction with the radius vector) to the stratification in the direction in which 
n increases. That is. 


Vt{U^{p) ; 0 <;. < 1 ( 2 . 1 ) 

and 

U;{V.(-H-) :0</i<l (2.2) 

represents the specific intensity of the ray travelling in the direction p and U~ represents 
the specific intensity of the ray travelling in the opposite direction. A finite set of values of 
: 1 < J < yri; 0 < ^1 < ^2 < ^3 • • - < 1) is selected 

or 


/ 



Unipi) 

Un{p2) 


\ 


\ UniPm) ) 


(2.3) 


and 
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U"(n) 



Figure 2.1; Schematic diagram showing a spherical shell with incident and emergent intensities; 
E~ or denotes contribution from internal sources. 
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(2.4) 


are m - dimensional vectors on Euclidean space. 

• The incident intensity vectors are (7^ and U~^i. 

• The emergent intensity vectors are U~ and Ui^i. 

The emergent radiation field will have the contributions from the internal sources say, 

E+(n + l,n) and E“(n, r» + 1) corresponding to the output intensity vectors U^+i and U~ 
respectively. 

Certain linear operators are assumed which reflect and transmit the incident radiation namely, 
f(n + l,n)f(n, n + l),r(n, n + 1) and r(n + l,n). These operators are calculated based on the 
physics of the medium- Then the output intensities are written in terms of the transmitted 
and reflected input intensities together with the internal sources as 


U~ = r(n + 1, n)Un + f(n, n 4-1), Un+i + ^ (”-7 ” + 1) (^•^) 

The introduction of the internal source terms namely,E+(n 4-1, n) and S“(n, n + 1) is due to 
Grant and Hunt [19]- The relationship given by equations 2.5 and 2.6 is called the Interaction 
Principle. This can also be written concisely as 


where 



5(n,n4- 



ui 

f^n41 



+ 1 ) 


5(n,n 4-1) = 


t(n + l,n) r(n,n + l)\ 
r(n4-l,n) t(n,n4-l)/ 


(2.7) 


( 2 . 8 ) 


2.2.2 Star Products 

If there is another shell with boundaries (n4-l}n4-2) adjacent to (n,n4"l), interaction principle 
for this shell can be written as 
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(2.9) 


( [;|» ] = S(n + l,n + 2) ( ^1; ) + E(n +1, n + 2) 

where S{n -|- l,Ti + 2) is similarly defined, as in equation 2.8. If one combines the two shells 
(n,n+l) and (n-l-l,n+2) then the interaction principle for the combined shell is written as,(for 
the thickness is arbitrarily defined) 

) =S'(n,n + 2)( ) +S(n + l,n + 2) (2.10) 

Redheffer [61] calls 5(n, n + 2) the star product of the two S-matrices 5(u, n + 1) and S(n 4- 
1, n 4 2); and 5(n, n 4 2) written as 

5'(n,n4 2) = 5(n,n4 l)**S'(n4 1,« 4 2) (2.11) 

Equation 2.11 is obtained by elimination and from equations 2.9 and 2.10. The r 
and t operators for the composite cell are written as 

t(n-l-2,n) = <(n + 2,n4-l)[r-K«+2.«4l)r(n,n+l)]"^i(n4 l,n) (2-12) 

t(n,n4 2) = i(fi,r» 4 1)[1 - r(n, n 4 l)r(n 42, n4 1)]"^<(»» 4 l,n + 2) (2.13) 

r(n + 2,n) = r(n4 l,») 4 t(n,n 4 + 2." 4l)r(n,n4 l)]“V(n4 2,n+1) (2.14) 

r(n,n4 2) = r(n 4 1, n 4 2) 4 <(« 4 2, n 4 1)[/- r(n,n+l)r(n 4 2, n 4 1)]“M«." +1) 


and 


S(n, n 4 2) = A(n, n 4 1; n 4 2)S(n, u 4 1) 4 A'(n; n 4 1,« 4 2)E(n 4 1, n + 2) 


(2.16) 


where I the identity matrix and 


A(n,n4l;n42) = 


A'(n;n4l,n42) = 


t(n 4 2, n 4 !)[/’-r(n,n 4 l)»*(n 42, u 4 1)]“^ 0 

t(n, n 4 l)r(n 4 2, n 4 !)[/ - r(n, n 4 l)r(n 4 2, n 4 1)] I 


I f(n 4 2, n 4 l)r(n, n 4 !)[/ - r(n 4 2, n 4 l)r(n, n 41)] ^ 
0 f(n,n4l)[/-Kn4 2,n4l)r(n,n4l)] ^ 
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and 


= ( 2 . 19 ) 

Similaxly, + 1,» + 2) is defined. 

To obtain physical interpretation of 2.12 the operator inverse is expanded in a power series. 
For example, 


OO 


f(n + 2 ,n) 

= SM« + 2,n) 

fc=0 

( 2 . 20 ) 

4 (n-l- 2 ,n) 

= t(n + 2,n4-l)[r(n,n-l-l) X 

( 2 . 21 ) 


r(n -h 2, n 4- l)]*f(n + 1, n) 

( 2 . 22 ) 


This operator acts on intensities to the right and gives the contribution to from U^. 
The term, tk{n + 2 ,n) may he recognized as diffuse transmission from n to n + 1, diffuse 
reflection from the layer (n, n+ 1 ), k times in succession and finally diffuse transmission through 
(n 4 -1, u + 2). Thus t{n + 2, n) is the sum of contributions involving scattering of all orders 
fc = 0,1,2,,.. OO. A similar interpretation can be given for other operators. 

E 5(a) is written to designate the cell, then 

5(a*/?) = 5(a)*5(/3) (2.23) 

where omr /I denotes the region obtained by putting the two cells a and 0 together. If the cells 
are homogeneous and plane parallel then 


(2.24) 

In general, the star multiplication is non-commutative. However, star multiplication is asso¬ 
ciative. For adding several layers a, /?, 7 ... 


5[(o(^ (/? ★ 7 ) Tfc*...)] — S[(^otk /3)k 7 ■*•...] (2.25) 

E the medium is homogeneous and very thick then one can use what is known as doubling 
method (see van de Hulst [30]). For example, 

5(2^rf) = 5(2'’-* d) * S(2'’-‘d), (P = 1 , 2 , 3 ,...) (2.26) 
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wldch means that we can generate the S-matrix for a layer of thickness 2^d in P cycles starting 
with S(d) rather than in 2^ cycles of adding the S(d)’8 one by one. If P = 10, then only a 
fraction 10/2^® ~ 10“^ of the computational work is needed to add 2^® layers of thickness d. 

2.2.3 Calculation of Radiation Field at Internal Points 

One expects the reflection and transmission operators to be non-negative on the physical 
grounds that intensities are always non-negative. This condition will be satisfied only when 
the optical thickness of the shell is less than certain value called the critial size or To.it- If 
the optical thickness r of the sheE in question is larger than the then we can divide the 
shed into several subshells whose r is less than the Tcrit a-ad then star algorithm is used to 
calculate combined response from the subshells whose total thickness is T, If, for example one 
need the radiation field at internal points in the atmosphere, one will have to divide the entire 
medium into as many shells as needed and c^culate the radiation field at the N points in 
the medium. One can write down the interaction principle for each shell and solve the whole 
system of equations (see Grant and Hunt [18]). 

Figure 2.2 shows the atmosphere in which to calculate the internal radiation field. The atmo¬ 
sphere is divided into N shells (homogeneous or inhomogeneous) with A and B as the inner 
and outer radii. The solution aJid £/“ (for any shell between shell l(at B) and shell N 
at A are obtained from the relations 


+ (2-27) 

tc = t(T.,n + l)t;-+i + V-n + l/2 (2.28) 

with the boundary conditions = U~{A). 

The quantities r(l,n -f 1), and are calculated by employing the initial conditions 
r(l, 1) = 0 and The computation is done by the following recursive relation 


r(l,nH-l) = r(n, n-f 1) + f(n ■+l,n)r(l,n)[/— r(n -I- l,n)rl,n]“^t(n,n-h 1) (2.29) 
W = r(n + l.u)V;ti/,+T„H/ai:-(n,n + l) (2-31) 


where 


t(n + l,n) = f(n +l,n)[/ —r(l,n)r(nH-l,n)]~^ (2.32) 
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and 


(2.33) 

(2.34) 

(2.35) 


r(n + l,n) = r(n +l,n)[/-r(l,n)r(n +l,n)]“' 
/in+i /3 = t(» +l,n)r(l,n) 

^'n+i/a = [I- r{n + 1,n)r(l, n)]"^ 


i(n, n + 1) = Tn+i/ 3 t(n, n + 1) (2.36) 

To calculate the radiation field at the internal points one proceeds as follows: 

1. Divide the medium into a number of shells (say N) with N-\-l boundaries as mentioned 
earlier. 

2. Start calculating the two pairs of reflection and transmission operators 
r(n4-l,n),r(n,n + l),f(n+l,n) and i(n,n4-l) in each shell (if the optical thickness 
in each shell is larger than Tcru then apply star algorithm to use doubling procedure 
if the medium is homogeneous). 

3. With the boundary condition that r(l,l) = 0 and = u‘''(6) and the r and t 
operators mentioned in subsection 2.2.2, compute recursively r(l,u -f- l),V ^/2 and 
t(n, n + 1) given in equations 2.29 to 2.30 from shell 1 to shell N (i.e.from B to A in 
figure 2.2). 

4. Next one sweeps back from A to B (see figure 2.2) calculating the radiation field given 
in equations 2.27 & 2.28 with the boundary condition U~^i = U~(a). 

One has to retain the operators r(l,n + l),t(n, n + l),V^j ^2 Ki+i/a she// that 

are calculated in the steps 2 and 3 above until one staHs calculating the radiation field from 
A. (The storage of these operators will not increase the memory of the machine required for 
computations. The operators can be stored on a magnetic disc and can he recalled whenever 
necessary. The operators can be removed from the disc or tape so that the disc or tape can be 
used to store fresh data). 

If the surface at A is reflecting, we can write 

^iv+-i = (2.37) 

where ra is the reflection operator. For a totally reflecting surface r© = I. Therefore, we have 


CGf=i = [r-'-(l,«' + l)rol-'V'+„^, (2.38) 

from which one can calculate from equation 2.37. Rest of the calculations follow equations 
2.27 suid 2.28. 
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2.3 Solution of the line transfer equation 


’I’liis section deals with the method to calculate the spectral line profiles using the Discrete 
Space Theory technique. The description given here follows from Peraiah [52, 53, 54]. 

In order to describe the method to solve the line transfer equation with partial frequency 
iredistribution functions and velocity fields, a two-level line transfer problem in non-LTE case 
in an extended spherical medium is considered. The equation of line transfer in such medium 
is given by 


gj(x,/i, r) dl{x,fi,r) 

^ dr r dn 


Ki\l3 + fi, r)][5(x, li, r) - I{x, fx, r)] 


and for the oppositely directed beam. 


(2.39) 


dl{x,-tx,r) l-fx^ dl{x,-ix,r) 
^ dr r dfi 


Ki\l3 + 4>{x, -IX, r)][5(x, -fx, r) - I{x, -fx, r)] (2.40) 


■where /(x,/i, r) is the specific intensity at an angle cos”^^(/i G (Ojl)) ^.t the radial point r 
and frequency (x = (i/ — i/q)/^*? where A, is some standard frequency interval). /? is the ratio 
KcJki of opacity due to continuous absorption per unit interval of x to that in the line. The 
source function S{x,±.n,r) is given by. 


^(x, ±/i,r) = 


</>(x, ±ix, r)S'i(x, r) -4- fdScjr) 
13 -I- 4>{x, ±ix, r) 


(2.41) 


where 5/(x, ±^,r) and Sc{r) refer to the source function in the line and in the continuum 
respectively. The line source function is given by 


Si{x, ±ix,r)= ^- T f dx' f R{x, ±fx] x\ ix')I{x\ /)d/ + eB{r) (2.42) 

^ ' 4>{x,±ix,r) J-oo J-i 


where e is the probability per scatter that a photon will be destroyed by collisional de-excitation; 
Sc 3 Jid B axe the continuum source function and Planck function, respectively; and (f>{x, r) 
and R[x, ±/i; x',/i') are the profile and partial redistribution functions. The profile function (f> 
is given by 


47r y y R{u, ft] I/', = <^(i/) (2.43) 

liquations 2.39 and 2.40 are solved by using the source functions defined in equations 2.41 
and 2.42 together with the appropriate redistribution function and following the procedure of 
integration described in Peraiah & Grant [57] and Grant h Peraiah [21]. This gives 
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Afm [U^n+l - + -^m^tTn+l/a] + ^n+l/^C/^ + ” 

Tn+l/3(/>)? + f^«)n+l/3^n+l/3* + |ni+l/3fl^n+l/2^i.n41/3^f'=-l“‘,n+l/3(**)^^^ 

Similarly, for the opposite-directed beam, 

Mm{U[^n-k-l - ^iT") ■•■ ^=[^^«7n+l/3‘^ + "^n+l/sC^ + ^i)n+l/2C^i,„+.l/2 = 

'TnJtXllipP + e<^i)n+l/35j,+i/3/» -|- ^ni+i/2Cr„+i/2^i.n+l/2S[.__iai,n+l/2(/l^)^^ 

where 

i 

^.>+1 = ^vrrJ^.1 1 (xi, -h/im, rn+i), 

^m,i,n+l /2 ~ +/im,’', 1 + 1 / 2 ), 

the A’s are the curvature matrices; pc is the curvature factor; and 


Bn+i/2 — 

C = [C>n.<^m/e]) 

A/m — [/^m^mA:]? 


(2.46) 


in which the /x’s and c’s are the roots and weights of the angle quadrature. The quantities <j) , 
etc. are to be similarly understood. The index n -f 1/2 refers the average over the shell 
with boundaries r„ tuid rn.+i, and i and i' are the indices corresponding to the frequencies of 
the incident and scattered photons. 

We can define coefficients WT;n+i /2 as 

i<f>iWk) = fli-,n+l/2Ci (2*47) 


where the subscript k is 

(i,i) =k = j^{i-l)JA<k<K = IJ 
with I and J being the total number of frequency and angle points, respectively. 
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Let 


K+1J2 = [^fcfc']«+l /2 = (/5 + <it)n+l/2Skk>, 

•®n+l/2 ~ ^<^fc)n+l/2^4.1/2^fcfe'5 

where T denotes transpose, and rewrite equations 2-44 and 2.45 as 

- P+] +p 4A+I7+^,,, +A-y;V,/,] + r„«/,*„+i/,t7+,,„= 
r„+l/J*»+l/J + J + V*(R) + P“ + C'"l„+,;3 (2 A8) 

and 

-t^n'+l/J -Pc[^'^Un+Ui + ^~^n+lj2] + r„+i/ 2 fli»+l/ 2 t^;Vl/j = 

T„+l/jJn+l /2 + \ [T(r<fKl>'^W] , 14 . 1 / 2 [t/'*' + U'^ +U~ + U~] „ 4 .i /2 (2.49) 

The average intensities U~^^j 2 iii the above equations are replaced by tbe following relations 
(Grant and Hunt [18]) 


{I — X^-^-Xl^Wn +-^n+l/ai^Ti+l = ^n+l/J (2.50) 

(/ — XnJ^\l2)U~ + Xn+l/al^n = ^n+1/2 (2.51) 

with X -\I for diamond scheme. Then by comparing the resulting set of equations with 
those in the canonical form of Peraiah & Grant [57] 


\ f <(n + l,n) 

t^n ) {r(n + l,n) 


r(n,n + l) \ ^ ( 
f(n,n+l) ) V ^^n+i 


+ 


^n+l/3 

^n 41/2 


) 


(2.52) 


we obtain two pairs of reflection and transmission operators, namely, t(u + lj^)» r(n,n -f-1), 
r(n + l,n) and t(n,n + 1). The expressions for these operators are given below. 


The cell operators are: 

t(n+ l,n) 
t(n,n+ 1) 
r(n + l,n) 
r(n,n+ 1) 

and the cell source vectors are: 


G-+[>l-i? + flr-V“], 
G^-g^-[T-\-A-D] 


^Ul/2 = A-S; 4 .i/ 2 ]t„+ 1/2 

^n+l/2 “ ^ ^ft+1/2 d" 9 


(2.53) 
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where 

G+- = [/-s+-s-+!-‘, 

G-+ = [/-s'^S+-]-S 
9+- = ir.+v>A+V'., 

g-* = iT„«/,A-n., 

B = Af —-r^+i/ 3 ^-, 

A = Af — 

A+ = [M + ir„+l/2Z+]"^ 

A“ = [M + —Ti^+i/aZ-l 

Kj. = (pA'^/r„+i/2) +i<y(/r‘*‘W^"'')n+i/2, 

y_ = -(M'/'J^n-fl/a) + ^^(^"^'‘’''Wi/25 
(T = 1 — C 

and I is the unitary matrix. 

These operators describe the radiation field in a moving medium and, therefore, all four redis¬ 
tribution functions 


R{x,-\-ii]x',n') (= 

Rlx,+fi]x',-ii') (= ), 

Rlx,-H]x\ii‘) (= ir+), 

Rlx,--H]x%-fi') (=ii"") 

are not equal and have to be calculated at each radial point of the medium. However, if 
the medium is static and partial redistribution is considered, then, by the symmetry relations 

(Hummer [31]), we have 

i 2 (a;,+/i;x', V) = -/i'); 

jR(x,+/i; x',-/i') = R{x^—H]x\-\-fAy,R:^~ = R'^ 

Furthermore, if one considers a static medium with complete redistribution, then all fom 
redistribution functions that evaluate the scattering integral reduce to approximately <f>(f> , 
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where <f> is tke profQe function and T is the transpose of the vector. Consequently, the form 
of the quantities Z+,Z_,l+,y_,etc., given above will be reduced to those given in Grant k 
Peraiah [21]. 

To obtain a stable solution one must choose an optical depth Tn+ 1/2 in each cell, given by 


rn+i/2 < T-crrt = mirifc 


iXk hp^tk 


(2.54) 


for the diagonal elements of the matrices A'^ and A”, and for the oflF-diagonal elements we 
must have 





< mink 

Tntiifc=v-n 

# . ... 

^tk' 


(2.55) 


Condition 2.56 is difficult to satisfy, which imposes a severe restriction on the size of the 
curvature factor to be used in each cell to obatin non-negative r and t matrices. Prom the 
inequalities set out in 2.55 and 2.56 above, it is clear that the medium of interest should be 
divided into a number of shells and relations (5.1) to (5.4) of Peraiah k Grant [57] should be 
used to calculate the diffuse radiation field at any point in the medium. However, it is time 
consuming to use a large number of sheila and subdivide each shell so that the curvature factor 
p is small enough to facilitate the use of the doubling algorithm described in Peraiah [50] and 
the star algorithm of Grant k Hunt [19, 20]. 

If the shell is halved p times, the star algorithm is repeated p times, and in this event the 
curvature factor and the optical depth for the subshell are given in terms of p, and r, 
of the shell 


» p.2-''/[l - p.(2-'-2-ni 

r„ = 

r' = + + + (2-56) 

where corresponds to a subshell approximately midway in the shell and p, is the curvature 
factor for the whole shell, defined as 


p. = Ar/r.., (2-57) 

being derived on the assumption that the optical depth in the whole shell is uniforni. is 
taken to be the mean value for all subshells, which introduces error, and to reduce this error 
one must divide the 3 hell into finer subshells that are small enough not to introduce serious 
errors. T is the mean radius of the subshell; R is the outer radius of the shell in terms of the 
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inner radius of the medium; and K — 2“' — 2~*. The relations set out in 5.57 are derived on 
the basis of equation (30) of Grant [17]. 

To conserve flux, one must ensure that the scattering integral given in equation 2.42, and the 
corresponding discrete equivalent given in the last brackets of equations 2.44, 2.45, 2.48 and 
2.49 are calculated exactly. For this purpose, the redistribution functions , etc., must 

be normalized to the machine accuracy so that, when there is pure scattering in the medium, 
conservation of flux can be checked easily (see Peraiah & Grant [57]). This can be achieved 
through the identity 

+ = l (2-58) 

^ P=l Q=l 


where 


Wp,Q = aiCj,a{ = 


AjRpQ 


and 


(^,Q) = i + (i+ tV 

and also by relation (4.3) of Peraiah &: Grant [57]. 

The method described above has been used to solve the line transfer equation in 
the following chapters. 
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Chapter 3 

Effects of angle-averaged Ey and R/j 
on emergent flux profiles in a static 
spherically symmetric stellar 
atmosphere 

3.1 Introduction 


The influence of partial frequency redistribution (PRD) is expected to become more pronounced 
in the extended atmospheres of supergiant stars because of low densities. Therefore, in this 
chapter, the dependence of the emergent flux profiles on the form of the angle-averaged partial 
redistribution functions Rv and Rjj is examined for spherically symmetric stellar atmospheres 
in a comparative way 

Scattering by a resonance line, that is broadened by radiation damping is described by the 
partial frequency redistribution function Rjj derived by Hummer [31]. This redistribution 
function is strongly coherent in the line wings leading to the lowering of line profile outside the 
Doppler core (Hummer [32], Vardavas [73]). Vardavas [71, 72, 73] studied the effects of angle- 
averaged and angle-dependent J?;, Rjj and Rm on spectral line formation in planar geometry 
for both static and moving media. Peraiah([51, 52]) studied the effects of partial frequency 
redistribution on spectral line formation in both static and expanding spherically symmetric 
stellar atmospheres using angle-averaged JR/ and Ru. Rangarajan et al. [60] studied the effects 
of stimulated emission on radiative transfer with PRD functions Rn and Rm in plane-parallel 
geometry. Milkey fe Mihalas [41] used a combination of JR// and Rm to study Solar Lyman-a 
line profile. 

Redistribution of photons in angle and frequency due to scattering in subordinate lines is given 
by Rv{x\n*\x,n) in the notation of Heinzel. Here x' and x represent the frequency displace- 


24 



ment from Kne center (in units of Doppler width) of the incident and the scattered photons 
respectively; n* and n give the directions of the incident and the scattered photons respectively. 
Heinzel [23] derived the correct form of angle-dependent laboratory frame redistribution func¬ 
tion (LFR), Rvj for radiatively broadened upper and lower levels of the atom. This LFR was 
extended to include coUisional broadening of the levels by Heinzel and Hubeny [24]. Hubeny 
and Hei n zel [29] solved the transfer problem using Ry for isothermal plane-parallel (both finite 
and semi-infinite) atmospheres. Mohan Rao et al. [44] examined the effects of Rn, Rm and 
Rv on source functions (5x) for plane-parallel and isothermal stellar atmospheres and showed, 
that, for a purely scattering medium with frequency independent incident radiation, the fre¬ 
quency dependent source function Sl{Rv) lies below Sl{Riii) but above *Sj[,(J?n) in fhe line 
wings. 

All the works cited above (except Peraiah [51, 52, 55, 56]) axe limited to planar geometry. The 
assumption of a plane-parallel atmosphere holds good only when the density scale-height in 
the atmosphere is small compared to the radius of the star. But, many stars, in particular the 
supergiants are known to have extended atmospheres; in the first approximation it is assumed 
that such stars are spherically symmetric. These stars are known to have low gravities leading 
to low densities in their extended outer layers where the PRD effects are expected to to be more 
pronounced (Basri [9]). In this study it is shown, the effects of Ry on emergent flux profiles 
are shown for the first time in comparison to the effects ofRjj and Ry considering spherically 
symmetric stellar atmospheres. The appreciation of these effects will be of relevance to the 
problems of spectral line modeling in stellar atmospheres and also to the modellers of stellar 
atmospheres. 


3.2 Equations and computational procedure 

The line transfer equation for a two-level atom in spherical symmetry is written as follows 
(Peraiah [52], Singh [65]) 


dl{x,±,,, r) ^ 1 aiix ±ix,r) ^ 4^x)][S{x, r) - /{i, ±1^, r)] (3.1) 

or r op 


where ± stands for the oppositely directed beams of radiation, /(x, ±/i, r) represents specific 
intensity of the ray making an angle cos~^fj, {p € [0,1]) with the radius vector, fl = 
(continuum to line opacity). 4>{x) represents the profile function, given as 


^(x) = / Riiy{x',x)dx' (3.2) 

OO 

where Rii{x',x) and Rv{x',x) are the angle-averaged partial frequency redistribution func¬ 
tions. Here, x' and x represent the frequency displacement from line center in units of Doppler 
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width of the incident and scattered photons respectively. The profile is normalized as 

/ +«» 

dx^(x) = 1 

-CX) 


(3.3) 


In the transfer equation (equation 3.1), S{xy r) is the total (Hne plus continuum) source function 

which is written as \a/ ^ , / 3 C / \ 

_ 4>(z)S,(x,r) + 0S,{r) ,3 . 

<Kx]+l3 

where Sc is the continuum source function set equal to 1 in this computation. Si{x,r) repre¬ 
sents the line source function, which , for a two-level atomic model, is given by the following 


expression 


5'/(a:, r) = ^TT f x)J{x')dx' -f- cB(r) 

<p[x) J~oo 


(3.5) 


where J{x') denotes the frequency-dependent mean intensity and eis the probability per scatter 
that a photon is destroyed by collisional de-excitation written as 

<^21 /O fl\ 


C21 + A2i[ 1 - erp(-^)]-i 


(3.6) 


where C-a. is the rate of collisional de-excitation from level 2 to 1 , A 21 is Einstein’s spontaneous 
emission probability for transition 2 —► 1, h is the Planck’s constant, u is the photon frequency, 
k is the Boltzmann’s constant and 0 is the temperature. 

Under the assumption of complete redistribution, the line source function is written as 

5,{r) = (1 - e) <l>{xf)J(x')dx' + £B(r) (3.7) 

J —CX) 

where <f>{x^ is taken to be a Voigt profile with damping parameter a. 

The line source functions written above for PRD (equation 3.5) and for CRD (equation 3.7) are 
valid for a two-level atomic model; whereas the simplest model needed to study subordinate 
line formation is a three-level atomic model- The reason for using a simplified two-level atomic 
model in this study is that, since we are interested in a comparison between radiative trans¬ 
fer solutions with three redistribution functions rather than in solving any specific transfer 
problem, we neglect the transition 3 1 (which provides an extra channel for destroying 

and creating photons in the subordinate line 3 2). Although this is a rather unphysical 

assumption that amounts to neglecting the terms A 31 , C 31 and Qe (rate of elastic collisions), 
it is suitable for our purposes. Under this assumption the line source function for a three-level 
atomic model reduces to a line source function for a two-level atomic model. 


3.2.1 Angle-averaged redistribution functions 

In this study, an idealized stellar atmosphere is considered assuming it to be isothermal, static, 
spherically sy mm etric and extended. The outer layers of this extended atmosphere will he 
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tenuous, where the scattering processes will be important. Also, the radiation field will be 
anisotropic and one should ideally use angle-dependent form of redistribution functions. How¬ 
ever, in the present study, the angle-independent (or angle-averaged) redistribution functions 
have been used. Although, by doing so, one loses information about the angular effects on 
emergent flux profiles, still, one is fully accounting for the frequency redistribution effects, 
which affect crucially the photon escape-probability, and hence the thermalization process. 
Therefore, by using angle-averaged redistribution fimctions, one accounts for critical aspects 
of redistribution process and sacrifices information only in an area of secondary importance. 

The angle-averaged JRjj has the following form (Mihalas [38]) 

Ru{x', x) = r e-“" [tan-\^) - (3.8) 

where x = max(|xl, |x'|) and x = mm(jx|, |x'l) ; a is the damping parameter. 

Rii envisions an atom having a sharp lower level and a radiatively broadened upper level. 
The above expression has been derived assuming isotropic phase function. Figure 3.1 gives 
the plot of the emission probability (i2j/(x',x)/<^(x^)) at frequency x per absorption, when the 
absorption is at frequency x\ Damping parameter a = 2.0 X 10~®. 

It is seen that a photon incident in the line wing has a high probability of being re-emitted 
in the line wing. Thus R[j is strongly coherent outside the line core. This property has 
implication on the line wing emissions from strong resonance lines. 

The angle-dependent Ry is given by the following expression (Heinzel [23]) 


Rv{x'^r^]x,n] = 


1 rtr/ 9 x-\-x' 9^ 

$ X — x' 6 


X Hiaicsc-,— 


csc-)-f Evix',x,9)] 


(3.9) 


where 6 is the angle between the incident (n’) and the scattered (ri) photon directions, a,- and 
aj are the damping parameters of the lower and upper levels of an atom respectively. The 
function Ey is given by 


where 


^ sin(9f2) ^ _x3 , _ 

Ey(x^x,9)- ^ R£ e (e 

.2wt 

(3.10) 

9 9 

A(t) = D{z -f icos-T -f aisec—) 

It £» 

0 

— D[z -{-tcos-) 

Jj 

(3.11) 

.x-fx\ 9 
z = (aj t 2 


(3.12) 

to = a,- -H aj — ix 


(3.13) 
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'I / /'. I 


Figure 3.1: Plot of the emission probability frequency x per absorption, 

when the absorption is at frequency x'. Damping parameter a = 2.0 X 10“®. D.Mihalas, Stellar 
Atmospheres, Chapter 13, p429. 





w' 

— Ci -f aj — ix' 

(3.14) 

D(w) 

= H{p,q)-{-iK{p,q) 

(3.15) 

UJ 

= p-iq 

(3.16) 


Here H{p,q) and K(jp,q) axe the Voigt functions which axe computed using the method of 
Matta & Reichel [37]. The angle-averaged expression for Rv is obtained by 

Rv{x\x) = 8it^ [ Rv{x',x^6)3in9d0 (3-17) 

This function is computed using a highly accurate computer program developed by Mohan Rao 
et al. [44] based on the method described by Heinzel fe Hubeny [25]. It is derived for an atom 
having lower and upper levels broadened by radiative broadening. The emission probability 
J?v(x', x) redistribution function has been plotted in figure 3.2 with damping 

parameters a* = Oj = 10"®. 

It is seen that a photon absorbed in the line wing has a high probability of being emitted at 
the line centre as well as in the wing. This is unlike the emission probability of Ru where, if 
a photon is absorbed in the line wing, it has the least probability of being emitted at the line 
centre. Also, Rv is less coherent than Rji in the line wing. 


3.2.2 Choice of depth points N 

The stellar envelope is divided into N depth points of equal radial thickness. These N depth 
points bound N — 1 spherical shells, each of equal radial thickness, n is the depth point counter 
that is taken to be increasing inwards and 


1 < n < iV (3.18) 

The outermost shell is bounded by n = 1 and n = 2; the innermost shell is bounded by 
n = A — 1 and n = N. T is the total optical depth of the atmosphere taken to be 10“*. 
We set iV = 75 in these computations. In selecting the value of JV = 75, the dependence 
of the emergent flux profiles on this parameter has been investigated to ensure the accuracy 
of results. In order to find out which value of N is suitable for this study, test runs were 
made with N = 10,25,50,75 and 100 and the emergent flux profiles were compared. Figure 
3.3 shows the emergent flux profiles due to CRD, Rn and Ry for various values of N using 
exponential opzw:ity distribution. 

Figure 3.4 shows similar plots for power law r"* opacity distribution. 

These plots show the sensitivity of the emergent flux profiles to the number of depth points 
N. It is seen that there is a lot of difference between the results due to iV = 10 and N = 25. 
This difference decreases for higher iNT = 50 and finally vanishes for iV = 75 and 100. This 
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Figure 3.2: Plot of the emission probability ■> ^)I ^'t frequency x per absorption, 

when the absorption is at frequency x'. Damping parameter a = 2.0 x 10 See P.Heinzel, 
1981, JQSRT, 25, 483. 


1 







X 


Figure 3.3: Sensitivity of emergent flux profiles to tte number of depth, points N using expo¬ 
nential opacity distribution. 
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Figure 3.4: Sensitivity of emergent flux profiles to the number of depth points N using r 
opacity distribution. 
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implies tkat, the resulting profiles axe quite sensitive to the value of N < 50, and insensitive 
to the change in N beyond N — 50. Therefore, for the parameters used in this study, N = 75 
is a suitable value. 


3.2.3 Optical depth laws 

For spherical atmospheres of fixed total optical thickness T, the optical depth variation (T(r)) 
in the atmosphere at radial points r is allowed to vary using 

1. exponential opacity distribution, and 

2. power law opacity distribution. 

Both these laws are appropriate for spherical geometries: the former being more realistic for 
static media; latter for constant velocity winds. Let us define to be the radius of the 
stellar photosphere. Further, let = 1 to bp the unit of length scale. Then, the extent of 
the atmosphere is defined in units of = 1. The atmospheric extent has been denoted by 
a variable R. For example, J? = 1 leads to plane-parallel approximation and il = 50 will 
mean that the atmosphere extends as far as 50 times the photospheric radius of the star. For 
plane-parallel cases the optical depth per layer is taken to be equal to be 2 T/N (see Kunasz 
&: Hummer [36]). For extended atmosphere the following laws have been used. 

From the definition of optical depth, we have 

’■(r.) = f''Kr)dr (3.19) 

Jrn. 

where k{r) is the opacity and is radius ( in units of 3U*) at the depth point n. 


Exponential opacity distribution 

For an exponential opacity distribution, following expression is adopted in this study 


A:(r) = (3.20) 

In this equation, q is a constant which controls the steepness of opacity variation. It is taken 
to be 0.1 in this study; ko is a constant which is evaluated by integrating the equation 3.20 
over whole radial length i.e. between ri and and setting the left-hand-side equal to T. It 
is foimd to be 


K = - 1 )-^ 


( 3 . 21 ) 
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Using the above relations, the optical depth of a shell bounded by n and »+1 has the following 
form 

) - T(r,,) = - 1)"^ (3.22) 

Power law opacity distribution 

For a power law opacity distribution we have (see Kiinawz & Hummer [36]) 

k{r) = kcr~"^ (3.23) 

where m = 0,2, 3. In the present study we set m = 2. ko is constant evaduated by integrating 
equation 3.20 (as done for exponential opacity distribution) over the whole radial extent of the 
atmosphere. It written as (with m = 2) 

ko = T(1 - (3.24) 

The optical depth of a shell bounded by n and n -f 1 is given by the following expression 

7-(r„„) - r(r„) = - r^') (3.25) 

It is seen that for exponential opacity variation, the matter is far well spread out in the 
atmosphere, whereas for opacity distribution where the matter is concentrated in the 
innermost shells. 

Following are the boundary conditions used for the transfer equation 

U+(x,T = 0,/i) = 0 and = T,fi) = Q if e > 0 (3.26) 

U+(x,T = 0,/z) = 0 and U;^+i(x, r = r,/i) = 1 if e=0 (3.27) 

wtidrc 

= 47rr^/(x,^,±T(r„)) (3.28) 

I{x, /t,T(r„)) being the specific intensity. + specifies a ray directed towards the bottom of the 
envelope (r = T), and - specifies a ray directed out of the envelope. The Planck’s function has 
been set equal to 1 throughout the medium. 

The above equations have been solved using the discrete space theory method 
described in chapter 2 of this thesis. 
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Figure 3.5: Optical depth variation versus depth points in a spherically symmetric atmosphere 
having = 50. n = 75 denotes the innermost boundary and n = 1 denotes the boundary of 
the outermost layer of the atmosphere. 






3.3 Results 


The results are presented for two sets of calculations. In the first set, the exponential opacity 
distribution given above is used and in the second set l/r^ variation of opacity given by 
equation 3.23 (by setting m = 2) is used. In each of these sets there are two subsets of 
calculations. In the first subset, the damping parameters for computing Rv are a,- = aj = 10“^ 
for the lower and upper level of the atom respectively, and in the second subset a,- = 10"'* and 
aj = 10"®. The corresponding damping parameter for computing Rn and CRD is taken to be 
a = o, + aj. Thus, for the first subset a = 2 x 10"® and for the second set a = 1.1 x 10"®. For 
convenience, the first set is called as ’model A’ and the second set as ’model B’. Also, the first 
subset is named as ’model a’ and the second subset as ’model b’. Using this nomenclature, the 
combinations are called as models ’Aa, Ab, Ba or Bb’. 

Several test runs were made in order to ensure the accuracy of this code besides checking the 
dependence of the accuracy of results on N as described earlier (figures 3.3 & 3.4) This code 
is able to reproduce the same line profiles as given by Kunasz & Hummer [36] in figure 8 of 
their paper and those given by Hubeny & Heinzel [29] in figure 2 of their paper. Moreover, 
flux conservation is maintained to an accuracy of the order of 10 in double precision for a 
purely scattering medium i.e. when e = 0 (Peraiah [56]). 

First, the results using exponential opacity distribution are presented. Figure 3.6 gives the line 
source functions of ’model Aa’. 

The upper set of curves is for a plane-parallel atmosphere R = 1 and the lower set is for an 
extended atmosphere with R = 50. The line source functions Sx for the extended atmosphere 
lie much below those for a plane-parallel atmosphere. This is the effect of atmospheric extension 
which increases the effective escape probability of photons; and the major effect of increased 
photon escape is to decrease Sx. Prom this figure we see that for a plane-parallel atmosphere, 
the line source function due to CRD always lies below those due to PRD. This is not the case 
for extended atmosphere where Sx due to CRD lies above those due to PRD for x > 4. Also, 
it is seen that for R = 50, Sx due to Ry always lies between Sx due to CRD and Rn- It is 
not so for a planar atmosphere where Sx due to Ry lies above both Rji and CRD for x > 4.5. 

Figure 3.7 gives the emergent flux profiles corresponding to the source functions given in figure 

3.6. 

An important result clearly seen is that in transition from a plane-parallel to an extended 
atmosphere the effects of Rjj and Ry are significantly enhanced (note the log scale). It is also 
seen that the flux profile due to Rn is much more affected than that due to Ry. The emergent 
flux profiles for R = 50 follow the vaxiation of the corresponding line source functions of figure 

3.6, whereas it is not so for a plane-parallel atmosphere. The emergent flux profiles due to Rjj 
and Ry lie below those due to CRD because both Rji and Ry are coherent in the line wings. 
Emergent flux due to Rn lies even below that due Ry because in the line wings Rn is much 
more coherent than Ry (see figures 3-1 fe 3.2). 
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Figure 3.6: Emergent line source functions for the parameters shown in the figure. 
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Figure 3.7: Emergent flux profiles corresponding to the source functions of fig.3.6. 
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Figure 3.8 gives the results of ’model Ab’. In this case, the lower level of the atom has been 
made ten times sharper (for Rv)', from a value of a,- = 10~® to a,- = 10“"*. 

There is a considerable change in the Sx curves due to Rv as compared to the curves given 
in figure 3.6. The values of Sx due to Rv have moved closer to the values of Sx due to Rn 
but have moved away from those of CRD. Thus the differences between the results due to Rv 
and CRD have increased. This is the consequence of making the lower level sharper and thus 
moving towards a condition for a resonance line where a,- tends to zero. 

Figure 3.9 gives the emergent flux profiles corresponding to the source functions shown in figure 

3.8. 

First by looking at the plots for j? = 1, it is seen that an increase in the ratio aj/cii from 1 in 
figure 3.7 to 10 in this case, has increased the departures between the emergent flux profiles 
due to Rv and CRD. This result is in qualitative aggrement with the result of Hubeny 
Heinzel [29]. In the line wings, the profiles due to Rv are closer to the profiles due to Rn 
than they were in figure 3.7. As mentioned above, this is the consequence of making the lower 
level sharper. Another noteworthy point is that these departures are enhanced by going from 
plane-parallel case (i? = 1) to an extended spherically symmetric envelope having R = 50. 

Now, the results using power law r"^ opacity are presented. 

Figure 3.10 gives the results for ’model Ba’. Significant differences between the results of this 
figure and the corresponding results of figure 3.6 are not found, except that for R = 50 case 
the Sx curve due to Rv lies closer to that of R// just outside the line core. Also in this figure 
R = 50 curves are much flatter than the corresponding ones in figure 3.6. This must be the 
effect of changing the opacity structure from exponential to r~^ form. 

Figure 3.11 gives the emergent flux profiles corresponding to the source functions in figure 
3.11. By comparing this figure with figure 3.7, the difference in the profile shapes for R = 50 
are clearly seen. This is the effect of changing the opacity structure. But more noteworthy 
is the point that in the far line wing the departures between the PRD and CRD profiles are 
the same as for exponential opacity distribution. Therefore we can say that in the far wings 
of the emergent flux profiles, whereas the atmospheric extent is responsible for the departures 
between the solutions due to PRD and CRD; opacity structure is not. It only changes the line 
shape for the kind of atmospheres considered here. 

Figure 3.12 gives the results of ’model Bb’. As compared to the results of figure 3.10, the 
line source function due to Rv lies closer to the Sx due to Rj/. As in figure 3.8, this is again 
the consequence of making the lower level of the atom sharper by a fax:tor of 10. As seen in 
figure 3.10, the curves for R = 50 are flatter than the corresponding curves (figure 3.8) using 
exponential opacity distribution. 

Figure 3.13 gives the emergent flux profiles corresponding to the source functions given in 
figure 3.12. The effect of Rv is seen to be significantly enhanced as was the case in figure 

3.9. Although the shape of the emergent profiles for R = 50 is conspicuously different from 
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Figure 3.8: Same aa fig. 4a but for flj/a,- = 10. 
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Figure 3.11: Emergent flux profiles corresponding to the line source functions of fig.3.10. 


43 





Figure 3.12: Same as fig.3.10 but for ajjai 







Figure 3.13: Emergent flux profiles corresponding to the line source functions of fig.3.12. 
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Table 3.1: Flux ratios showing the departures between the emergent fluxes due to PRD and 
CRD at a: = 4, for different values of atmospheric extension R. R — I represents plane-parallel 
approximation, whereas R>1 represents an extended and spherically symmetric atmosphere. 
These results are for aj/ai = 1,10 



aj/ai = 1 

aj/ai = 10 

R 

F^{CRD) 

_ 

K(UHD) 

FJRv) 

F4CRD) 

FsiRrr) 

FARv) 

1 

d.l54E+00 

0.221E-f01 

0.122E-h01 

0.120E-I-00 

0.261E-1-01 

0.215E-I-01 

5 

0.517E-01 

0.351E-1-01 

0.141E-1-01 

0.339E-01 

0.407E-f01 

0.306E-f01 

10 

0.290E-01 

0.397E-f01 

0.146E-f01 

0.185E-01 

0.456E-1-01 

0.331E-f01 

50 

0.691E-02 

0.488E-1-01 

0.153E-h01 

0.422E-02 

0.547E-I-01 

0.375E-f01 

100 

0.440E-02 

0.500E-1-01 

0.154E-h01 

0.267E-02 

0.558E-1-01 

0.381E+01 


Table 3.2: Same as above but for a: = 6 



aj/ai = 1 ‘ 

Oj/ai = 10 

R 

F^iCRD) 

FACRD) 

KICHU} 

~my) 

F^CRB) 

F,IURD] 

TM-rA 

Kiumn 

1 

b.987E-01 

0.146E-h02 

0.152E-1-01 

0.674E-01 

0.242E-f-02 

0.623E-f01 

5 

0.249E-01 

0.510E-f02 

0.173E-f01 

0.153E-01 

0.730E-t-02 

0.847E-h01 

10 

0.134E-01 

0.681E-1-02 

0.176E-H01 

0.817E-02 

0.921E-1-02 

0.879E+01 

50 

0.298E-01 

0.104E-f02 

0.178E-I-01 

0.178E-02 

0.126E-f-03 

0.915E-i-01 

100 

0.187E-02 

0.112E+02 

0.179E+01 

0.112E-02 

0.132E+03 

0.918E-I-01 


the shape of corresponding profiles given in figure 3.9, the differences in the far wing between 
the PRD and CRD profiles do not seem to be different. Thus again, the change in opacity 
structure does not change the departures between PRD and CRD in far wings. 

In order to further elucidate the above effects, the folio wing, three tables are presented. These 
tables give the ratios of the emergent fluxes due to CRD to the corresponding emergent fluxes 
due to PRD {Rj[ and Rv), indicating the departures of the results of PRD from those of 
CRD at n = 4,6 and 8 for different values of the atmospheric extension R. Along with these 
ratios, CRD results are explicitly tabulated so that the ratios can be unravelled. Because the 
emergent flux profiles due to PRD and CRD are nearly the same in the line core (r G [0,3]), 
and most of the departures are in the line wings, we chose to present the emergent flux ratios 
only for a; = 4,6 and 8 in these tables. As defined earlier, iJ = 1 denotes the plane-parallel 
approximation, whereas, R> I denotes extended spherically symmetric atmospheres. 

Tables 3.1, 3.2 and 3.3 gives the results at a: = 4, ® = 6 and at x = 8, respectively. The 
comparison is shown for the ratios aj/ui = 1 & ay/a,- = 10. The effect of Rn is clearly more 
than the effect of Ry ,both as a function of x and R. The given flux ratios are seen to increase 
much more sharply in going from R = 1 to 5 than in going from R = 5 to higher values, 
especially for Rjj cases tabulated here for x = 6 and x = 8. Thus the geometry of the system 
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Table 3.3: Same as above but for x = 8 



aj/oi = 1 

R 

F*(CRD) 

FJTCHET 

F^lRii) 

f/THn 1 

1 

0.653E-01 

0.149E-I-03 

0.166E-1-01 

5 

0.144E-01 

0.706E-H03 

0.179E-1-01 

10 

0.764E-02 

0.815E-I-03 

0.180E-1-01 

50 

0.166E-02 

0.789E-f-03 

0.181E+01 

100 

0.104E-02 

0.781E+03 

0.180E-1-01 


F^(CRD) 


dilai — 10 

“TT^- 


rmu] 

rn 




0.417E-01 

0.868E-02 

0.458E-02 

0.988E-03 

0.602E-03 


0.330E+03 

0.957E+03 

O.lOOE+04 

0.876E+03 

0.860E-}"03 


0.857E+O1 

0.960E+01 

0.968E-}"01 

0.973E-HO1 

0.972E-f01 


plays an important role for resonance line transfer in tbe far wings (x > 6). At frequencies 
X = 4 and 6 the departures between PRD and CRD results keep on increasing from = 1 to 
100, but at X = 8 the maximum departure occurs at R = 10 for R/j and at R = 50 for Ry for 
the physical parameters considered in this study. It is obvious from the results for the ratio 
flj/a,- = 10, that the departures between PREf and CRD results are considerably more than 
those in the case where aj f Oi = 1. Increase in R enhances these departures in the same trend 
as when Cj/ai = 1. 


3.4 Conclusions 


The effects of angle-averaged Ry in spherically symmetric stellar atmospheres are investigated 
in comparison to the effects of angle-averaged Rjj and complete redistribution approxima¬ 
tion. This study establishes that for static and isothermal atmospheres, the partial frequency 
redistribution effects are significantly enhanced in transition from a plane-parallel to an ex¬ 
tended spherically symmetric stellar atmosphere. In a comparative way, it has been shown 
that sphericity affects the solutions of R/j more than the solutions of Rv- The effect of Ry 
on the emergent flux profiles increases substantially by making the lower level of the atom 
sharper. This effect is found to get further enhanced in transition from a plane-parallel to an 
extended atmosphere cause. These differences between the emergent flux profile due to PRD 
and CRD do not seem to vary by altering the opacity distribution in the atmosphere, at least 
for the exponential and forms used in this study. In the line wings, the departures between 
the results due to PRD and CRD increase much more rapidly for small atmospheric extensions 
than for larger extensions. 
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Chapter 4 

Partial frequency redistribution effects 
of Ry and Rjj in expanding spherically 
symmetric stellar atmospheres 

t 


4.1 Introduction 


In the last chapter, we presented the comparative partial frequency redistribution effects of Ry 
and Rii on the spectral line formation in spherically symmetric static stellar atmospheres. In 
this chapter, we examine the same problem by including a constant velocity of expansion. 
In the previous chapter, the velocity fields were neglected in order to estimate the effect of 
sphericity alone. This chapter shows the combined effect of sphericity and velocity field, 
using the partial frequency redistribution functions Ry and i?//, in comparison to the effects of 
the assumption of complete redistribution (CRD). Thus this study inspite of being parametric, 
goes a step closer to a realistic situation of extended atmosphere and velocity of expansion. 
The atmosphere in assumed to be inhomogeneous and isothermal. 

Mohan Rao et al. [45] examined the effects of Rj/, Rm and Ry on the line source functions 
using planar geometry and static isothermal atmosphere. They showed that, for a purely 
scattering medium with frequency independent incident radiation, the frequency dependent 
source function Sl{Rv), Hes below Si(Rin) but above Sl(Rii) in the line wings. Singh [66] 
investigated the comparative effects Ry, Rn and CRD on the emergent flux profiles, using a 
spherically symmetric and static stellar atmosphere. Hubeny k Heinzel [29] solved the transfer 
problem using Ry and Ru for isothermal finite and semi-infinite atmospheres. Besides these 
few works, it is difficult to find a comparative study in literature involving Rv, Rn and CRD. 
Therefore, in this chapter, we present, for the first time, a comparative study of the partial 
frequency redistribution functions Ry, Rn and CRD showing the combined effect of sphericity 
and velocity fields. 
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As higliliglited in the previous chapters, this study is an interesting study of astrophysical 
relevance, because the PRD effects are expected to he pronounced in the extended atmospheres 
with low surface gravities. The physical parameters used in this study are characteristic of the 
atmospheres of cool giants and supergiant stars. 


4.2 Equations and computational procedure 


H and x are the frequency displacements from line center, in units of standard Doppler 
width, of the incident and scattered photons respectively, seen in the rest frame, then, the 
corresponding frequency displacements, seen in the fluid frame at radius r, are 


X' = x'±V{r)fi (4-1) 

X = x ±‘V{r)pi (4*2) 

where [0,1]) is the cosine of the angle between the radius vector and the direction of 
propagation of the radiation, db stands for the oppositely directed beams of radiation, and 
y(r) gives the macroscopic velocity at radius r. In this chapter we assume that ^the 
atmosphere expands with a uniform velocity, i.e., we neglect velocity gradient 
across the atmosphere. This kind of situation may arise in circumstellar shells where one 
can neglect the velocity gradients by assuming a constant wind velocity. 

The equation of line transfer in spherical symmetry, rest frame and for a two-level atomic 
model is written as (Peraiah [52], Singh [65]) 

^ l^ a/(x ^ ^ ^)] (4.3) 

dr r Ofi 

where i stands for the oppositely directed beams of radiation, /(a:, i/ij r) represents specific 
intensity of the ray making an angle cos~^fi € [0,1]) with radius vector at the radial point 
r. (/>[x, ±/i, r) represents the profile function, given as 

<l>{x,±ii,r) = <f){Xyr) = I Rij,v{X',X)dX' (4-4) 

J-OO 

where angle-averaged paurtial frequency redistribution functions RnyiX , X) have been stated 
and explained in the previous chapter (see equations 3.8 & 3.9 and figures 3.1 & 3.2). 

Redistribution function Rv is derived for subordinate line scattering and assumes an atom 
having both the upper level and the lower level broadened by radiative and collisional damping. 
Rii envisions an atom having a sharp lower level and a radiatively broadened upper level. 
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S{xj ±/x, r) is the total (line plus continuum) somce function which is written as 


S(x, ±fi,r) 


(^(X.r)5|(a:.±/i.r)+/35c(r) 

<f>{X,r) + p 


(4.5) 


where (3 = Kc/Ki (continuum to line opacity, Sc is the continuum source function set equal to 
1 in this computation. 5(r, i/i, r) represents the line source function, which is given by the 
following expression 


S,{x,±^,r) = j*^d^'R„,v(X',X)I(x',±i,,r) + cB(r) (4.6) 

where c is the probability per scatter that a photon is destroyed by collisional de-excitation. 
The procedure to compute c for the chosen resonance lines has been given above in section 5.3. 

f 

Under the assumption of complete redistribution, the line source function is written as 

5^(r) = l^y^ dx dfi<i>{X,r)I(x,±iJi,r)-{■ eB{r) (4.7) 

where (f>{X, r) is taken to be a Voigt profile with damping parameter set equal to the sum of 
the damping parameters (for Ru and CRD) of the upper and lower levels of an atom causing 
subordinate line scattering. 

The normalization condition for the profile function at each radial point is 

/ -foo 

dx<f>{X,r) =: 1 (4.8) 

-OO 

In a differentially expanding medium, due to the presence of velocity gradients, on has to 
compute, at all radial points all the four redistribution functions appearing in the scattering 
integral viz 


Rii,v{x' -t- V"(r)/i, X 4- V{r)n) 

Rii.vix' - V{r)(j, X - V{r)fji) 

Rii,v{x' — V(r)/i, X -+■ V(r)/i) 

in order to evaluate the diffuse radiation field (see chapter 2 &: Peraiah [52, 55]). 
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Figure 4.1: Emergent flux from an expanding atmospkere having the velocity of expansion 
V = 1 in units of thermzd velocity; i? = 1 denotes plane-parallel case and R — ^Q denotes the 
extended atmosphere having an extension 30 times the stellar radius (set equal to 1). 
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Figure 4.2: The effect of expansion velocity V = 1 in comparison to static case V — 0 for an 
extended atmosphere R — 30. 
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In this work since we assume a constant expansion velocity, we do not have velocity 
gradients across the atmosphere. This assumption is valid for expanding circumstellar shells; 
and reduces the computational time of the Jiv" and Rji matrices enormously. The above written 
redistribution functions are computed at one depth point in the atmosphere, and the other 
corresponding redistribution functions at other depth points are equal to these. As is known 
from the works of Bemat [11] Sanner [63] and several other authors, the ratio of expansion 
velocity to the turbulent velocity for cool giants and supergiants lies between 0.5 and 2.0. 
Therefore, in this study for all the cases the value of this ratio has been chosen to be equal to 
1.0 at all depth points of the atmosphere. The envelope is divided into n shells, n = 1 defines 
the outermost shell and n = N defines the innermost shell. The total optical depth,T, of the 
medium has been set equal to 10'*. The total number of shells have been taken equal to 25 
m this study. The sensitivity of the results for static cases has been discussed in detail in the 
previous chapter section 3.2.2. Because at this stage, we are interested only in the departures 
between the solutions due to various redistribution functions, rather than modeling any specific 
spectral line, the choice of iV = 25 is justified. 

♦ 

Ideally, for moving media calculations, one should use angle-dependent redistribution functions, 
because the frequency-angular coupling effects could be important in anisotropic radiation field. 
By employing angle-averaged redistribution functions, we have neglected the angular effects 
(which are of secondary importance to this study) but have fully accounted for frequency 
redistribution effects. Also, it has been shown by Vardavas [73] that for planar geometries and 
differentially expanding medium, the differences between the results due to angle-dependent 
and angle-independent redistribution functions are below 10%. Moreover we have considered 
only constant velocity of expansion leading to the absence of velocity gradients, therefore, it is 
justified to use angle-independent redistribution functions. 

Following boundary condition has been used for the transfer equation 


where 


U^{x,fjL) = iirrll{x,fi) (410) 

being the specific intensity, -f specifies a ray directed towards the bottom of the 
envelope (r = T), and - specifies a ray directed out of the envelope. 

To solve the line transfer equation the discrete space theory method described in chapter 
2 of this thesis has been used. 
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4.3 Results 


Figure 4.1 shows the differences between the PRD and CRD results for a plane-parallel and 
extended atmosphere expanding with a constant velocity at V = 1 (in terms of thermal ve¬ 
locity). The differences are substantial for both the cases. Closer look at the profiles 
shows that the differences are slightly more for the extended atmosphere. Recalling the results 
of the previous chapter, we know, that sphericity enhances the PRD effects. It is also known 
&om earlier works (e.g. see Hempe [26]) that increasing velocity decreases the PRD effects. In 
the present case, where both extension and velocity is present, the effects cancel out and the 
differences remain more or less the same for jR = 1 and for R = 30. 

Figure 4.2 shows the effect of introducing a constant wind velocity Y = 1 (in thermal-velocity 
units) into an extended atmosphere R = 30. Asymmetry in the profile shape introduces 
enhanced departures between corresponding PRD and CRD results just outside the line core 
X w ±3. Overall differences remain significant. 

Figure 4.3 shows the effect of increasing the extension from R = 30 to (over three times) 
R = 100. Differences are still found to be substantial. Another effect well known effect of 
increasing the extension is seen in this figure. The central absorption becomes shallower- This 
is because the effective emitting area increases in the line core by increasing the atmospheric 
extent. 


4.4 Conclusion 


The role of partial frequency redistribution functions Ry and Rji on the emergent flux profiles 
has been presented in a comparative way. Extended atmosphere and constant wind velocity 
has been assumed. It is shown by this study, that for the parameters used here, the effects 
of partial frequency redistribution are substantial for all cases discussed here. These lesu ts 
clearly demonstrate the importance of partial frequency redistribution functions in spectral 
line formation in the atmospheres of cool giant and supergiant stars. 
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Chapter 5 


Effects of angle-averaged Rju on 
spectral line formation in expanding 
spherically symmetric stellar 
atmospheres 


5.1 Introduction 


After having explored the combined effects of sphericity and velocity fields on the spectral line 
formation, using partial frequency redistribution functions Rv and iZj/, we now consider the 
redistribution function Rni alone. The motivation behind doing so is the following. Most 
often the limiting case of complete redistribution (CRD) is taken to be an adequate repre¬ 
sentation of the partial frequency redistribution (PRD) function Rni iii radiative transfer 
problems to study spectral line formation in stellar atmospheres. The works of Finn [16] an 
Vaxdavas [71, 72] have shown that for plane parallel (static and moving) stellar atmospheres 
the differences between the emergent intensities due to Rm and CRD are below 20 o. ut, 
with the advent of very high resolution stellar spectroscopy (resolving power )’ 

highly axicurate line profiles (e.g. for bright giants and supergiants see Banner [63] are aval 
able, therefore, even small differences between CRD and PRD profiles could be important in 
a quantitative analysis of stellar spectra e.g. in stellar winds and mass loss phenomenon in 
giants and supergiants (Hempe [26]). Also, the combined effect of sphericity and macroscopic 
velocity on these differences cannot be estimated a priori. 

So far, there is no paper in literature showing the effects of Rm alone on spectral line for¬ 
mation in spbericaJIy symmetric and differentiaJIy expanding stellar atmospheres. Therefore, 
tins paper has considered a parametrized spherically symmetric and differentially expanding 
stellar envelope to study the effects angle-averaged Rm on spectral line formation. Line trans¬ 
fer equation has been solved in the rest frame of the star assuming two-level atom model. 



This work we have also examines the behaviour of emergent mean intensity profiles (7*), the 
source function (5) and the emergent flux profiles (1^) as a function of atmospheric extension 
(6/fl =outer/inner stellar radius), maximum velocity of expansion (Vj) and the thermalization 
parameter (c). The results are in aggrement with the earlier known results of Kunasz and 
Hummer [36] and Mihalas [38]. The mean intensity profiles due to Rin{x\x) and CRD we 
compared. Here x* and x are the frequency displacements from the line centre in Doppler umts 
of the incident and scattered photons respectively. 


5.2 Basic equations and computational procedure 


If X* and X are the frequency displacements from line centre, in units of standard Doppler 
width, of the incident and scattered photons respectively, seen in the rest frame, then, the 
corresponding frequency displacements, seen in the fluid frame at radius r, are 

X' - x' ±V(r)/i (5-1) 

X = a:±y(r)/i (5*2) 

where /^t(€ [0, Ij) is the cosine of the angle between the radius vector and the direction of 
propagation of the radiation, dz stands for the oppositely directed beams of radiation, and 
l^(r) gives the macroscopic velocity at radius r. 

The line transfer equation for spherical geometry in the rest frame for a two-level atom model 
has the following form (Peraiah citeperlh) 


dr r dfji 

where zt stands for the oppositely directed beams of radiation, /(r,zb/i,r) represents specific 
intensity of the ray making an angle cos~^fi (/i € [0,1]) with radius vector at the radial point 
r. represents the profile function, given as 

<j>{x,±jJL,r) = 4>{X,r) = I 

J —cx> 


where angle-averaged partial frequency redistribution function Riii{X'^X) has the following 
form (Mihalas [38]) 

X —— ) - tan~^{ ^ ~ -)]du (5-5) 

where a is the dzunping parameter set equal to 10“^ in this calculation. S{x, dzfi, r) is the total 
(line plus continuum) source function which is written as 


S{x,±^i,r) = 


(p{X, r)Si{x, jzp, r) -i- l3Sc{r) 


(5.6) 
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where = KcJKi (continuum to line edacity, Se is the continuum source function set equal to 
1 in this computation. 5(a:, represents the line source function, which is given by the 

following expression 

1 1 — C /+O0 /•f‘1 

= 2^7]^^ J_^ d^'Ri„{X',X)l(x',±i>,r) + €B{r) (6.7) 

where e is the probability per scatter that a photon is destroyed by coUisional de-excitation 
written as 

Q 

Ctx -f- A2 i[ 1 - exp(^)]-i ^ ^ 

where C%\ is the rate of collisional de-excitation from level 2 to 1, A 21 is Einstein’s spontaneous 
emission probability for transition 2 —>■ 1, h is the Planck’s constant, v is the photon frequency, 
k is the Boltzmann’s constant and 6 is the temperature. 

Under the assumption of complete redistribution, the line source function is written as 

1 — f /Hhoo /Hhl i 

Si(r) = ^ dx d,i^(X, r)I(x, ±/,. r) + eB(r) (5-9) 

where 0(X, r) is taken to be a Voigt profile with damping parameter set equal to 10“^. 


The normalization condition for the profile function at each radial point is 

r+oo 


/ +00 

dx<^(X,r) = 1 

-OO 


(5.10) 


In a differentially expanding medium, due to the presence of velocity gradients, on has to 
compute, at all radial points all the four redistribution functions appearing in the scattering 
integral (equation 5.9) viz 
R(x'-i-V(r)/i,x + V(r)/i), 

R{x' - V(r)p, X - V(r)p), 

R{x'-{■ V{r)iJ,,x — V{r)n) and 
R{x' - K(r)/i, x -f V(r)/i) 

in order to evaluate the diffuse radiation field (Peraiah [52, 55]) 

In this work, a linear velocity law given in Peraiah [52] has been used, which gives a radially 
increasing velocity, with minimum velocity (K) at the innermost boundary and maximum 
velocity of expansion (H) at the outer boundary of the stellar envelope. The envelope is 
divided into n shells, each of equal radial thickness, n = 1 defines the outermost shell (r = 0) 
and n ~ N defines the innermost shell (r = T). Total optical depth (T) of the envelope 
has been taken to be 10® and we have set = 10 in these computations. Optical depth 
(t) through the envelope varies as r“®. 14 has been set equal to 1 and 2 in units of mean 
thermal velocity of the gas. Rest frame calculation is limited to low expansion velocities due 
to numerical difficulties. Following are the boundary conditions used for the transfer equation 

U^{x,T = 0,/i) = 0 and U^^^{x, r = r,/i) = 0 if c > 0 (5.11) 

(7/'(x,t = 0,/i) = 0 and r = T,fi) = I if e = 0 (5.12) 


58 




T 


Figure 5.1i Source function S versus optical depth t through spherically symmetric and diifer- 
eatially expanding envelope having b/a = 5, 50; Vb = 0, 1; and T = 10^. Corresponding CRD 
results are shown by dotted (....) lines. 


wtlC/ITC 

= 47rr5[/(a;,^,r(rn)) 

I{x,n,r{rn)) being the specific intensity. + specifies a ray directed towards the bottorn of the 
envelope (r = T), and - specifies a ray directed out of the envelope. The Planck’s function has 
been set equal to 1 throughout the medium. Flux conservation has been maintained accurate 
to the order of 10”^® in double precision for purely scattering medium i.e. when c 0 (Peraiah 
[56]). 

To solve the above equations the discrete space theory method described in chapter 2 of this 
thesis has been used. 
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Figure 5.2: Emergent mean intensities in the line versus corresponding to source functions 
of figure 1. X = (i/ — v'(i)/w ; w being standard Doppler width. 

5.3 Results 

Figure 5.1 gives the run of the source function S' as a function of optical depth t through a 
spherically symmetric and differentially expanding envelope where e = /? = 10~®. We recover 
the known results that the increase in velocity and extension suppresses the source functions. 
We also see that the effect of expansion velocity is more at larger extension. The source 
functions due to CRD are leas than those due to Rni. 

Figure 5.2 gives the mean intensity profiles corresponding to the source functions shown 
in figure 1. We see that at any given velocity of expansion the mean intensity in the line is 
lower for larger extension. The profiles are asymmetric for Vi = 1 (expressed in mean thermal 
velocity of the gas). More important is the result that the difference between the CRD and 
Riii profiles is below 5% near the line centre and negligible at the wings. This is much smaller 
than reported for planar geometries by Finn [16] and Vardavas [71, 72], 
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Figure 5.3: Emergent line flux Fx versus x at Vj = 1; c — 10 6/a — 1 (plane-parallel 

atmospkere), 5, 50 and T = 10^. 
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Figure 5.4: Emergent mean intensity profiles versus x for a spherically symmetric envelope 
having extension hja — 50 ; expansion velocity Vh—Q (static case) and 1; and c = ^ = 10“^. 

In figure 5.3 emergent flux profiles aje shown for a differentially expanding (14 = 1) envelope 
at different values of extension. Emergent line flux from a plane-parallel (i>/a = 1) atmosphere 
is considerably higher than from spherically symmetric extended (6/a = 5, 10) atmospheres. 
Again, the known result is recovered that higher the extension, lower is the emergent line flux. 

In figure 5.4 a comparison is made between the emergent mean intensity (7*) profiles due 
to CRD and Rm for a static (14 = 0) and differentially expanding (14 = 1) spherically 
symmetric envelope having extension 6/a = 50 and e = /3 = 10”*. Although, the differences 
between CRD and Rju profiles are slightly more than what they were in fig. 2 where a lower 
value of e = = 10”® was used, still, the differences are within 5%, again, much lower than 

those reported for planar geometries by Finn [16] and Vardavas [71, 72). 

Figure 5.5 shows the effects of e and on emergent flux (F*) profiles due to Rm for a static 
and differentially expanding spherically symmetric and extended envelope. It is found that the 
profiles have deeper absorption depths and more asymmetry at lower value of e and (3. The 
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same is also true for CRD as skown by Kunasz and Hummer [36]. 


5.4 Conclusion 

Tke differences between the emergent mean intensity profiles due to complete redistribution 
aad angle-averaged Rm are foimd to be much smaller in spherically symmetric and differen¬ 
tially expanding stellar atmospheres than those reported for plane-parallel (static and moving) 
atmospheres. 
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Chapter 6 


Partial frequency redistribution effects 
on observed resonance lines using 

model atmospheres of cool stars 

» 

6.1 Introduction 


After considering the parameterized atmospheres of cool stars in the previous chapters, this 
chapter employs realistic computed model atmospheres of late type giants and 
supergiants to study PRD effects on some of the observed spectral lines. The 
spectral lines chosen for this study are the three strong resonance lines Si II A1818, C IIA1335 
and Ba II A4554. Out of these, the first two spectral lines have been observed in the cool 
supergiant star a-Orionis, spectral type M2 la-Ib, with very high resolution from the Hubble 
Space Telescope (see Carpenter et al. [12]) using the Goddard High Resolution Spectrograph 
(GHRS). The observed high resolution profiles of Ba II A4554 have been presented by Sanner 
[63]. The motivation behind this study is that PRD effects are expected to become important in 
the extended tenuous atmospheres of cool giants and supergiants. This is an interesting study of 
astrophysical relevance because it shows the model dependence of the differences between the 
results due to PRD and CRD in observed spectral lines. The knowledge of these departures 
could be useful to the modellers of stellar atmospheres and for quantitative analysis of the 
stellar spectra. 

As mentioned in the review of the related work (section 1.4), several researchers have made 
attempts, with varying degree of success, to model the chromospheres of late-type stars based 
on the matching of observed and computed PRD profiles. Most of these works have neglected 
velocity fields and atmospheric extension. In the present study both the atmospheric extension 
of the chromospheres and velocity fields have been taken into account. This study aims 
at highlighting the differences between the solutions of PRD and CRD and the 
dependence of these differences on the adopted model atmospheres. 
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6,2 Adopted model atmospheres 

In this study two kinds of model atmospheres of cool giants and supergiants have been adopted 


• atmospheres without chromospheres; and 

• atmospheres with chromospheres. 

The atmospheres without chromospheres have been taken from Plez [58]; and to study PRD 
effects in atmospheres with chromospheres, the extended chromospheric model of a-Orionis 
(M2 la-Ib) given by Hartmann and Avrett [22] has been adopted. 

The model atmospheres without chromospheres are the spherical opacity sampling models for 
M-giants and supergiants (Plez [58]) covering the range 3000 K < Tc// < 4000 K, —0.5 < 
^1-5 for 1, 2 and 5 Af©. These assumte hydrostatic equilibrium, spherical symmetry, 
LTE and homogeneity. Out of these given models, the models having the following parameters 
have been adopted in this study: 

1. T^jf = mOK'Jogg = - 1.0; R//?© = 166.0; and logL = 3.804. 

(Te//, 5 , M,R and L are the effective temperature, surface gravity, mass, radius and 
luminosity, respectively). 

2. Te// = 3000R:; logg = 0; M/Mq = 1.0; R/Rq = 166.0; and logL = 3.304. 

The run of the gas pressure Pg, the electron pressure Pe and the temperature T for the first 
model has been plotted in figure 6.1, and the same parameters for the second model have been 
plotted in figure 6.2. Hereafter we will call them as modtll and modelB. These parameters 
have been plotted against the variable n which is the depth point counter in the atmosphere, 
n = 1 refers to the outer most boundary and n = 25 refers to the innermost boundary of the 
atmosphere. 

The main difference in the two adopted models is that modell is cooler {Ttff — 3000 K), 
than the model2 (T*// = 4000 K). This choice has been made to study the effect of different 
temperature structure on the differences between the CRD and PRD results. Modell and 
model2 do not have chromospheric temperature structure. 

The third model atmosphere, hereafter models, chosen for this study is the extended chromo¬ 
spheric model of the supergiant star o:-Orionis (popularly known as Betelgeuse, spectral^ fyp® 
M2 la-Ib, given by Hartmann & Avrett [22]. This model has been constructed by a desire to 
explain cool winds observed from a wide variety of gaints. The authors have used a schematic 
wave-driven wind theory to compute the density, temperature and velocity structure of the 
wind. The run of the hydrogen number density nn cm~^, the electron number density cm 
and the temperature T have been plotted versus z = r/fi* in figure 6.3. z = 30 refers to the 
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Figure 6.2: Model2: Variation of the gas pressure Pg, the electron pressure Pc and the tem 
perature T inside the adopted model atmosphere (Plez (1992)). 
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point 30 stellar radii away from z = 1 . In this model steUar radius, i 2 « = 1350i{Q and stellar 
mass, jkT* = lOAf 0 . 

hord^ to conapute the emergent flux profllea using the above model atmospheres, one requires 
to solve the equation of line transfer; and to solve the transfer equation one needs to compute 
the thermalization parameter e and the line center optical depth tq, at all the radial points r 
inside the atmosphere. In the next three sections, it is explained that how one proceeds to 
achieve that. 


6.3 Calculation of thermalization parameter e 


The thermalization parameter e is defined as the probability per scattering that a photon will 
be destroyed by collisional de-excitation. As seeq. in earlier chapters, the line source function 
has two terms, the first term with coefficient (1 — e) is the noncoherent scattering term, (see 
equation 3.7), followed by a thermal source term cBp. The line source function will depart 
from being a Planckian (B,,) if c ^ 1 . In such cases the assumption of local thermodynamic 
equilibrium (LTE) does not hold. In virtually all situations of astrophysical interest, c << 1 
in regions of line formation. Only very close to the photosphere, one finds that e -+ 1, and 
only in such regions the condition of LTE prevails. Thus, when e « 1, the noncoherent 
scattering dominates the line source function, and the partial frequency redistribution (PRD) 
effects should be important. The expression to calculate the thermalization parameter at radial 
point r is written as 


e 



C2i{r) 


C 2 i{r) 4- Aaifl - exp 



( 6 . 1 ) 


where C^i is the rate of collisional de-excitation from level 2 to 1 at each radial point r, A^i is 
Einstein’s spontaneous emission probability for transition 2 —+ 1 , h is the Planck’s constant, 
M(i is the line center frequency, k is the Boltzmann’s constant and T(r) is the temperature at 
ladius r. In this study, (721 (r) and T(r) are depth dependent. r(r) at each depth is taken 
fiom the respective model, and C 2 i(r) is computed fr 6 m the following expression given by Van 
Regemorter [70] 


C«(r) = 20.60AL/V.(r)v/WP (6-2) 

where N^[r) is the electron number density (cgs units) at radial point r, taken from the respec¬ 
tive model at each depth point, AE is the energy of the transition and k is the Boltzmann’s 
constant. The function is interpolated from table 2 of Van Regemorter [70]. 
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Figure 6.4; Variation of thermalization parameter e Versus depth point CTunter n for Ba II A 
4554, plotted for modell {T,jf = 3000 K) and for model2 (Z,/ = 4000 K). 




Figure 6.5: Variation of thermalizatioii parameter e versus the atmospheric extent z for C II 
A 1335 and Si II A 1818, plotted for models (a-Orionis). 
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t?** 

figure 6.4 shows the variation of thermalization parameter e versus depth point counter n for 
® ^ 4554, plotted for modell (Tc// = 3000 K) and for model2 (r*// = 4000 K).It is seen 

at c < .25 < 1 in the outer layers of modell and model2. In these regions the noncoherent 
scattering donunates the line source function.. 

Figure 6.5 shows the variation of thermalization parameter € versus the atmospheric extent z 
or C IIA 1335 and Si II A 1818, plotted for models (ot-Orionis). fYom this plot it is clear that 
orr>l, c«i. This is the effect of the chromospheric temperature structure and the large 
^elision of the atmosphere. For most of this atmosphere, the situation is highly non-LTE. 

m principle, as compared to modell and model2, the effects of non-coherent scattering will be 
larger. 


6.4 Calculation of line center optical depth tq 


another important ingredient required for solving the line transfer equation is the line center 
optical depth r©. The expression to compute the line center absorption coefficient at radial 
points r in the transition between levels i and j is written as 


xM = 


rWeC Auz>(r) 


(6.3) 


Consider a shell bounded by radii r and r + Ar. If we assume that Xkq have a steep 

variation through this shell of thickness Ar, then To(r) will be given as 


7 b(r) = X Ar 


(6.4) 


where 


me = 9.1094 x 10-=®^ 
c = 2.9979 X 10^°cm/s 
k = 1.3807 X \Q~^^erg/deg{K) 

axe the mass of an electron, speed of light in vacuum and the Boltzmann’s constant respectively; 
ATj is the number density of absorbers in level n computed by solving the Saha’s ionization 
equation and assuming solar abundances taken from Allen [1]; fij is the oscillator strength for 
the transition between i and j; uq is the line center frequency; Az/£)(r) is the Doppler width, 
which is written as 
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where is the mass of the particular atom. 

Using the expressions given above, the line center optical depth is calculated for Ba II A 4554, 
C II A 1335 and Si II A 1818. To study PRD effects in non-dhromospheric atmospheres given 
by modell and model2, Ba II A 4554 has been chosen. The line center optical depth structure 
throu^ the atmosphere of modell (Te// = 3000 K) and model2 (71.// = 4000 K) is plotted 
in figure 6.6. It is seen that r© for model2 is higher than it is for modell. This is because 
of higher temperatures there are more number of singly ionized barium (Ba II) atoms in the 
atmosphere of model2. 

Figure 6.7 shows the variation of line center optical depths C II A 1335 & Si II A 1818 versus 
the atmospheric extension z in the chromosphere atmosphere of a-Orionis (model3). It is seen 
that To of C II A 1335 remains significantly abbve the unit optical depth and drops to Tq = 1 
at large extension z 20. This structure will have a bearing on the emergent flux causing a 
deep absorption self-reversal. The tq varia.tion of Si II A 1818 exhibits a plateau with Tq ^ 10 
at « « 11 slowly falling to a value of tq « 1 at z w 18. 


6.5 The equation of line transfer 


With the thermalization parameter e and the line optical depth Tq computed, one goes ahead 
to solve the equation of line transfer. Since strong resonance line axe chosen for this study, the 
partial frequency redistribution function function Rii{x\x) given by equation 3.8 and plotted 
in figure 3.1 of this thesis, is appropriate and has been used. This redistribution applies to an 
atom having radiation damping in the upper state and coherence in the atom’s rest frame. The 
results due to Rii{x\x) are compared with the results due to complete redistribution (CRD) 
approximation. 

Here, x' and x axe the frequency displacements from line centre, in units of standard Doppler 
width, of the incident and scattered photons respectively, seen in the rest frame. Then, the 
corresponding frequency displacements, seen in the fluid frame at radius r, are 


X' = x'±V{r)p (6.6) 

X=x± V{r)n (6.7) 

where /i(€- [0,1]) is the cosine of the angle between the radius vector and the direction of 
propagation of the radiation, ± stands for the oppositely directed beams of radiation, and 
V (r) ^ves the macroscopic velocity at radius r. 




Figure 6.6: Variation <>f ^ iddnC= 

n in tlie atmosplieres of modell ( e// 


the depth point counter 

4000 K). 



Log[T-p(z)] 


a-Orionis M2 la-Ib 


C n X1335 


^Si II \1810 


Figure 6.7: Variation of line center optical depths of C II A 1335 & Si 11 A 1818 versus the 
atmospheric extension z in the chromospheric atmosphere of a-Orionis (model3). 






Tke equatioa of line transfer in spherical symmetry, rest frame and for a two-level atomi 
model is written as (Peraiah [52], Singh [65]) 


dlix,±ti,r) l-iA^dIix,±fA, r) r)][5(ar, r) - f(x, r)] (6.8) 

^ dr r dii 

where db stands for the oppositely directed beams of radiation, /(x, ±/i,r) 
intensity of the ray making an angle co 5 ”^/i {fi € [0,1]) with radius vec or a e ra 
r. ^(a:, i/i, r) represents the profile function, given as 


<^(x,±/i,r) = <f>{X,r) = Rii{X',X)dX' 


(6.9) 


where angle-averaged partial frequency redistribution function Rn{X ,.X^) has th " g 

form (Mihalas [38]) 

where = max{\xl \x'\) and x = mtn(|xl, \x'\) ; a is the damping parameter. 

Ru eavisio,is aa atom having a sharp lower level and a f thlt 

above expression has been derived assuming isotropic ph^e function. ^ 

gives thTplot of the emission probahiUty [R„[x',xmx')) " I>“ ^t^o^-tion, 

when the absorption is at frequency x'. Damping parameter a . X 

5(x, ±^,r) is the total (Hne plus continuum) source function which is written as 


S{x,±n,r) 


(f>(X, r)5/(x, db/i, r) + PScjr) 

^iX,r) ^ 


( 6 . 11 ) 


where ^ = KJK, (continuum to line opacity. S. is the continuum “““ ““ 

1 in this computaLn. 5(x,±p.r) represents the hne source function, which is given by 

following expression 




( 6 . 12 ) 


where e is the probahiUty per scatter that a photon is destroyed by colUsional de^exdtation. 
^hr;r:.::dme to Xute efor the aosen resonance lines has be«i given above in section 5.3. 

Under the assumption of complete redistribution, the line source function is written as 
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(6.13) 


1 —. £ f+OO /+1 

^'(o=y oo / j ±/'. *■)+«*('■) 

where ,^(Jf,r) is taken to be a Voigt profile with damping parameter set equal to 10~®. 
The normalization condition for the profile function at each radial point is 



(6.14) 


As e^lained in the previous chapters, in a differentially expanding medium, due to the presence 
of velocity gradients, on has to compute, at all radial points all the four redistribution functions 
appearing in the scattering integral (equation 5.12) viz 


Rii{x' + V(r)/i, X + V{r)fx) 

Rii{x' - V(r)^, X - V(r)fA) 

Rn{x' + V(r)^, X - V{r)pi) 

Rn{x' - V{r)fi, x-^V{r)fi) 

in order to evaluate the diffuse radiation field (see chapter 2 h Peraiah [52, 55]). 

In this work, a linear velocity law given in Peraiah [52] has been used, which gives a radially 
increasing velocity, with minimum velocity (Va) at the innermost boundary and maximum 
velocity of expansion (Vt) at the outer boundary of the stellar envelope. For the model3, one 
could have computed the ratio of expansion velocity to the turbulent velocity from the values 
given in the table 2 of Plez [58]. It has not been done as the values of expansion velocities given 
are larger by a factor of two to three. As is known from the works of Bernat [11] Banner [63] and 
several other authors, the ratio of expansion velocity to the turbulent velocity for cool giants 
and supergiants lies between 0.5 and 2.0. Therefore, in this study for all the models the value of 
this ratio has been chosen to be equal to 2.0 at the outermost layer and decreasing towards the 
bottom of the atmosphere according to the linear velocity law mentioned above. The envelope 
is divided into n shells, each of equal radial thickness, n = 1 defines the outermost shell and 
n = N defines the innermost shell. The total number of shells have been chosen according to 
the number of depth points given in the respective model atmospheres adopted for this study. 
The line center optical depth is computed as explained in section 5.4. Following boundary 
condition has been used for the transfer equation 


tfti(x,/i) = 0 and Uj^^^ix.fi) = = 1) (6.15) 

where 
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= 47rrJ/(ar,/^,r(rn)) (6.16) 

being the specific intensity. + specifies a ray directed towards the bottom of 
the envelope (r = T), and - specifies a ray directed out of the envelope. = 1) is the 

Planck function at line center frequency i/q and at temperature T at n = 1. 

To solve the line transfer equation the discrete space theory method described in chapter 
2 of this thesis has been used. 


6.6 Results 


In this section, the emergent fulx profiles due to partial frequency redistribution function 
Ril{x\x) are presented in comparison to those due to the assumption of complete redistribution 
(CRD), in order to show the departures between the two. 

First the results due to the formation of Ba II A 4554 in non- chromospheric models, modell 
and model2, are presented. Figure 6.8 shows the emergent flux profile of Ba II A 4554 using 
modell having Te// = 3000 K. The profiles due to PRD (Rn) aJid CRD hardly differ from one 
another, except for negligible departure in the far wing on the red side. 

Figure 6.9 shows the emergent flux profiles of Ba IIA 4554 using model2 having Te// = 4000 K. 
We see « 10% differences between the profiles due to PRD and CRD. As compared to figure 
6.8, the departures between PRD and CRD are more in this case because the star is hotter 
and that leads to more number of singly ionized Ba atoms. 

Figure 6.9 shows the emergent flux profiles of Ba IIA 4554 using model2 having Te// = 4000 K. 
We see 10% differences between the profiles due to PRD and CRD. As compared to figure 
6.8, the departures between PRD and CRD are more in this case because the star is hotter 
and that leads to more number of singly ionized Ba atoms. This can be seen from figure 6.6, 
where throughout the atmosphere, the line center optical depth of Ba II is more for the hotter 
star. Thus, for this line, hotter the star, more will be the differences, with the upper limit 
being decided by the ionization potential of barium. 

Figure 6.10 shows the emergent flux profiles of Si II A 1818 due to PRD and CRD in the chro¬ 
mospheric atmosphere of a-Orionis. These profile exhibit enormous difference in the wings 
between PRD amd CRD profiles. This is a collisionally dominated line and exhibits the typical 
absorption self-reversal. One can clearly see how the assumption of CRD can lead to signif¬ 
icantly different (and erroneous) results. The assymetry in peaks is due to the presence of 
velocity fields. One can also notice that, as expected in the line core, there are no differences 
between the results due to CRD and PRD. 

Figure 6.11 shows the emergent flux profiles of C II A 1335 from the chromospheric atmosphere 
of a-Orionis, due to PRD and CRD. One immediately notices the striking differences between 
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Figure 6.8: Emergent flux profiles of Ba II A 4554 for modell. 




X 


Figure 6-9: Emergent flux profiles of Ba II X 4554 for model2- 
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Figure 6.10: Emergent flux profile of Si II X 1818 from the chromospheric atmosphere of 
a-Orionis, modeI3, due to PRD and CRD, 
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the results due to PRD and CRD in the line wings. Also as compared to the Si II profiles, the 
absorption in the line core is much deeper. Again, one can see that how much lowering of the 
line profile in the wings of this strong resonance line takes place by employing partial frequency 
redistribution function Rn- This is because^this function is strongly coherent in the line wings. 
In the process of quantitative analysis of stellar spectra, these substantial differences between 
PRD and CRD rAn lead to incorrect temperatute structure and hence incorrect estimates of 
the quantities dependent on temperature. It can be emphasized, on the basis of these results, 
that for line modeling and for computing synthetic spectra of chromospheric atmospheres of 
cool giants and supergiant stars, the assumption of complete redistribution should be relaxed. 


6.7 Conclusion 


The study of the differences between CRD and PRD emergent flux profile of strong resonance 
fines Ba IIA 4554, Si II A 1818 and C II A 1335^ from the atmospheres of cool giant and super- 
giant staxs has been presented. Atmospheres with chromospheres and without chromospheres 
have been considered. The differences between PRD and CRD are found to be highly model 
dependent. For non-chromospheric type atmospheres, the overall differences between the PRD 
and CRD profiles of Ba II A 4554 are within 10% for the hotter star having T^ff = 4000 K, and 
negligible for cooler star having T,ff = 3000 K; For the chromospheric model of cool super¬ 
giant star a-Orionis (M2 la-Ib), there are substantial differences between the PRD and CRD 
results for Si II A 1818 & C II A 1335 emergent flux profiles. On the basis of these results, it 
is recommended that the assumption of CRD should be relaxed for quantitative analysis 
of spectra of cool gaints/supergiants having chromospheres. 
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Chapter 7 

Concluding Remarks 


7.1 In retrospect 

In the last six chapters we reported on the effects of partial frequency redistribution functioM 
on the emergent flux profiles against the background of relevant hterature. J 

high resolution spectral data from satellites Uke the Hubble Sp^ 
and many others is making tUs an important area of study for the ^cmate qu 
analysis of spectra. This thesis is presented as a modest contribution to the efforts in that 

direction. 

Besides "true absorption”, scattering of radiation is one of the important physical proces^ 

in the stellar atmospheres. In very extended and tenuous atmospheres of 

sapergiant stars, where the rates of colHsions are low, scattering of photons is the 

factor in transport of radiation. Early works on radiative transfer Msumed that “ 

dther strictly coherent or that the photons are completely redistnbut^ over ^e p A 

Neither of these Umits is achieved exactly in stellar atinosph^es and 

use partial frequency redistribution functions. Due to the difflcdties in tack mg 8 

integral employing exact partial frequency redistribution functions, the numerical soluto of 

the transfer eiation had to wait tiU modern times with their effective comput^ri t»k 

The progress of studies in this direction is now rapid with varying degr^ “ 

important to be aware of the departures of the emergent flux profiles due to partial frequency 

redistribution from those due to complete redistribution. 

This thesis is an attempt to locate the various scenarios in the atmosphere of cool ^ 

the effects of partial frequency redistribution functions Rji, Rm and Rv (as mmp 
CRD) could be important and to present these effects in a comparative manner using spheri^ 
geometry. The study using parameterized stellar atmosphere provide an insight mto the 
role of single or combined effects of physical variables like the atmospheric extension, ve oci y 
fields andL thermaUzation parameter. The study involving reahstic model atmosphere and 
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observed strong resonance lines is an important step further in investigating the role of partial 
frequency redistribution functions on spectral line formation in cool giants and supergiants. 


The present work 

The assumption of complete redistribution is expected to fail in extended stellar atmospheres 
where low gravities and hence low collision rates prevail. In this thesis, with a clear objective of 
bringing out the departures between the solutions due to PHD and CRD in various situations 
in the atmospheres of cool giants and supergiants, the understanding has been developed in 
stages, starting with spherically symmetric static atmospheres, then including the velocity fields 
and finally considering realistic temperature and density structures to explore PRD effects in 
observed resonance lines. This systematic approach has resulted in considerable success in 
generating useful information and enhancing the understanding of the PRD effects on spectral 
line formation in cool giant and supergiant stais. 

• Chapter one explains the need to investigate this problem and presents the layout of 
the thesis. 

• Chapter two, describes the solution of equation of radiative transfer in the frame-work 
of the discrete space theory of radiative transfer using Peraiah-Grant formalism. 

• In chapter three, the effects of angle-averaged Rv and Ru on emergent flux profiles 
in static spherically symmetric stellar envelopes have been examined. The equation of 
transfer for spherical symmetry in the rest frame is solved for a two-level atom. Initially, 
to bring out the effects of geometrical extension alone, velocity fields have been neglected. 
A constant value of the thermalization parameter has been assumed and the Planck’s 
function is set equal to unity throughout the atmosphere. The background continuum is 
neglected. This study establishes that sphericity enhances PRD effects. It has also shown 
that the effect of spherical extension is much more on the emergent flux profiles resulting 
from Ru than on those resulting from Ry when the damping parameters of the lower 
and upper levels of an atom are comparable. The effect of Ry (as compared to CRD) on 
emergent flux profiles increases by making the lower level sharper than the upper level. 
This effect is further enhanced by increasing the spherical extension. 

• In chapter four, a step further is taken to include velocity of expansion in order to explore 
the combined effects of sphericity and expansion on the emergent flux profiles resulting 
from Ru and Ry, as compared to CRD. An atmosphere having geometrical extension and 
small velocity of expansion is considered which is characteristic of the atmospheres of cool 
giants and supergiants. It is assumed that there is no backgroimd continuum; Planck’s 
function is set equal to unity and a constant value of the thermalization parameter is 
chosen at all depths in the atmosphere. This study shows that there are substantial 
differences between the emergent flux profiles resulting from PRD {Ru and Ry) and 
CRD. 
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• la chapter five, the PRD effects of Rm alone on the emergent line profiles in an expand¬ 
ing spherically symmetric stellar atmosphere are explored. This is because, the limiting 
case of complete redistribution is usually taken to be an adequate representation of Rui 
in radiative transfer problems to study the spectral line formation in stellar atmospheres. 
However, the combined effect of sphericity and macroscopic velocity cannot be estimated 
a priori. Only low velocity regime is considered which is appropriate for the atmospheres 
of cool giants and supergiants. This problem is treated in the rest &ame formalism. In this 
work it is established that the differences between the emergent mean intensity profiles 
resulting from CRD and Rni *ire much smaller in expanding spherically symmetric stellar 
atmospheres than reported previously for plane-parallel atmospheres. 

• Chapter six presents the study of the differences between CRD and PRD emergent flux 
profile of strong resonance lines Ba II A 4554, Si II A 1818 and C II A 1335, from the 
atmospheres of cool giant and supergiant stars. Atmospheres with chromospheres and 
without chromospheres have been considered. The differences between PRD and CRD 
are found to be highly model dependent.; In non-chromospheric type atmospheres, the 
overall differences between the PRD and CRD profiles of Ba II A 4554 are within 10% for 
the hotter star having Te/y = 4000 K, and negligible for cooler stax having Te// = 3000 
K. There are found to be substantial differences between the PRD and CRD results for 
Si II A 1818 & C II A 1335 profiles computed using the chromospheric model of cool 
supergiant star a-Orionis (M2 la-Ib). On the ba^is of these results, it is recommended 
that the assumption of CRD should be relaxed for quantitative analysis of spectra 
of cool giants/supergiants having chromospheres. 


...and its limitations 


Although the foregoing achievement is significant, still we must look at the limitations of the 
currently achieved success, particularly with a view to plan future fine of action. In the present 
phase of work, we were eager to explore the frequency redistribution effects on the spectral 
line formation in the atmospheres of cool giants and supergiants. Even within the scope 
of employing exact forms of the redistribution functions R//, Rm and Rv, small velocity 
fields, atmospheric extension, use of computed model atmospheres to examine PRD effects on 
observed resonance lines of Ba II A 4554, Si II A 1818 and C II A 1335, we restricted ourselves 
in many directions. For instance, we considered only angle-averaged redistribution functions 
in this study. By doing so, we accounted fully for the frequency redistribution effects but 
sacrificed information on mgular effects, which could be important when the radiation field is 
anisotropic. We assumed a two-level atomic model for all cases of study. .Thus we neglected 
extra channels of excitation or decay. But since only strong resonance lines are considered, we 
could afford to neglect this aspect. We did not take up the modeling of these observed lines, 
but aimed at highlighting the differences between the PRD and CRD results using realistic 
model atmospheres. We did not adopt the velocity structure given by the model atmospheres, 
but employed a linearly increasing velocity field. The value of this velocity field matches the 
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values derived from the observational data. Considering the uncertainties in the model values, 
it is justified to employ the velocity structure used in this study. 


7.2 In prospect 

Further possibilities for exploring the partial frequency redistribution effects on spectral line 
formation are many and most of them can be handled by the techniques developed and em¬ 
ployed in the present work. Some of these are listed as follows: 

1. The radiation field in the tenuous atmospheres can be anisotropic where the angular 
effects on the emergent profiles could be significant. This would require the use of angle- 
dependent partial frequency redistribution functions. 

2. Keeping in mind the availability of very/high resolution spectroscopic data for stellar 
atmospheres the assumption of an isolated line may fail. One has to consider the use of 
multilevel radiative transfer employing angle-dependent redistribution functions derived 
for multilevel atoms. 

3. Modeling of chromospheric lines must relax the assumption of complete redistribution, and 
exact angle-dependent forms of redistribution functions should be employed. Attention 
needs to be devoted to the study of the combined effects of partial frequency redistribution, 
velocity fields and the geometry along with other physical phenomena. 

Some of these questions have been partially resolved, but much work remains to be done before 
a completely satisfactory understanding of spectral line formation in the presence of partial 
redistribution effects is achieved. 


A journey of thousand miles begins with a single step. 

— Chinese proverb 
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